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The task-oriented groups considered here consist of a number of individuals, each having 
initially one piece of information which must be transmitted to all the others to complete the 
task. Interest is centered on the communication net which restricts the possible channels for 
messages. At every sending time each individual sends all the information he has acquired to 
one other individual; the major assumption here is that this recipient is chosen at random from 
the possibilities given by the communication net. The information state is defined as a matrix 
which shows where the initial information has spread. These matrices can be considered as the 
states of a Markov chain, and in this way the distribution of completion times for the task 
is obtained. Some special cases are worked out and generalizations are indicated. 

A proof is given of the formula for the shortest possible completion time in any net when a 
fixed number of messages is sent by each individual at each sending time. 


1. Introduction. One of the outcomes of Kurt Lewin’s suggestions to- 
ward the group dynamics approach to psychology has been experiments 
on cooperative problem solving by small groups of people. The type of 
experiments with which we are particularly concerned here have been 
reported by A. Bavelas (1950) and L. S. Christie, R.D. Luce, and J. Macy, 
Jr. (1952; this will be referred to as CLM). 

The basic question of how people behave when they must work to- 
gether to solve a problem or complete a task is, of course, tremendously 
complicated and in these experiments attention is directed mainly to one 


~ aspect: the effect of various restrictions on the possibilities for communica- 


tion between individuals. The attempt is made to minimize the complexi- 
ties due to such factors as individual differences in intelligence and per- 
sonality, differences in motivation, etc. by making the task of trivial 
simplicity so that no ingenuity is required and there can be little emotional 
involvement. A typical task is as follows. Each individual is given a set of 
colored marbles, only one color being common to all the group. The mem- 
bers of the group must exchange messages about their own colors and 
what they have learned about the colors of the others, until finally every- 


body knows the common color. To eliminate the effects of differences in 
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speed of perception and response, messages are sent only after everyone 
has indicated readiness to transmit; the transmissions then take place 
simultaneously, each individual sending to just one of his possible re- 
cipients. Each member of the group knows initially to whom he can send 
messages, but does not know to whom the others can send. 

In the following we consider some mathematical questions arising from 
these experiments. One of the main aims is to determine what type of 
communication net (i.e., restrictions on the paths for messages) is opti- 
mum for reducing the number of messages required to complete the task. 
This question is not answered here, but it is hoped that the present con- 
siderations will provide a start toward the answer. Of course, even the 
greatly simplified simulation of the group problem-solving situation which 
the above experiments represent is still too complex for a complete mathe- 
matical treatment; important aspects must be omitted and further sim- 
plifying assumptions are needed. Thus the questions of group learning— 
to which much attention is devoted in CLM—and of morale and satisfac- 
tion are not included in this treatment. 

We refer to the number of messages that are sent by each individual to 
complete the task as the completion time, T. This amounts to taking the 
time between transmissions as the time unit. 

The major assumption here is that messages are sent at random, that is, 
at each sending time each individual picks a recipient from among those to 
whom he is permitted to send messages with equal probability* of picking 
any one of them. This is, of course, not a correct statement of how people 
will behave in small task-oriented groups, because it means that the same 
information might be repeated to a recipient, thereby failing to communi- 
cate new information, instead of deliberately choosing a recipient to whom 
the sender’s information may be new. In spite of the obvious objections to 
this assumption of randomness it seems worth while to consider this prob- 
lem for the following reasons. 

In the first place, although the assumption of random choice of recipient 
can not be exactly correct, it does seem to approximate the actual be- 
havior of the group at first, if the communication net is fairly complex. 
This is indicated by the results reported in CLM. The distribution of com- 
pletion times under the assumption of random sending (abbreviated to 
RD, random distribution, in the following) was computed approximately 
by a Monte Carlo method in connection with these experiments. It can be 
seen (chap. III of CLM) that for the more complicated communication 
structures the observed distribution of completion times for the first six 


* The assumption of equal probability is not essential as noted in Section 6. 
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trials was not too far from the RD. Of course, there was learning with 
increasing experience, so that later trials showed improvement over ran- 
domness. With simple communication structures the observed times de- 
parted considerably from the RD right from the start, but even in these 
cases the RD gave a rough indication of what was to be expected. 

The last remark points to a second use for the RD. This is as a basis for 
comparison, an indication of what to expect in any case. It is obviously of 
value to have such a basis even though divergence from it is expected. 
The divergence itself has significance as an indication of rational rather 
than random behavior. 

Finally the assumption of randomness can be justified on the principle 
that difficult problems can often be solved by treating simplified versions 
of the problem first and then adding the additional complicating factors 
needed to bring the situation closer to reality. The alternatives to the 
random assumption would require considerations about what rational 
behavior would be for each individual at each stage, and also to what ex- 
tent the actual individuals involved can be expected to be rational. These 
are obviously questions which are not easy to answer. The only relaxation 
of randomness which might not present too much difficulty is the mini- 
mum assumption that a sender will never repeat a message when he is cer- 
tain that it carries no new information to the recipient. 

2. Minimum completion time. We consider first a question not involving 
the assumption of random sending. For a group of given size, , and any 
possible arrangement of the communication net, what is the minimum 
completion time assuming now that the individuals are told where to send 
their messages so as to achieve this minimum. This minimum time, r, is 
the best that can be done under any conditions given m and the number 
of messages, 7, which each individual can send at each time. For the usual 
case of one message at each time, the value of 7 has often been stated but 
the writer does not know of any published proof. The minimum time, 7, is 
in this case the smallest integer containing log, . If we allow more than 
one, say 7, messages to be sent by each individual at each time, then 7 is 
the smallest integer containing log, 7, which is the same as stating that 7 
is the integer satisfying 
logrt41 M7 <1 +log,+i 1 \ 


n<(r+1)7*<(r+1)n. 


(1) 


or 


A proof of this result is given in the Appendix. For 7 = 1, the lower bound 


for + was established by Macy and Luce. 
3. Notation. From now on we consider the problem of Section Deuse 
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determination of the RD. To describe the structure of the communication 
net and the flow of information we use matrix notation (cf. Shimbel, 1951, 
and CLM). 
Let the individuals be numbered 1, 2, . . . 2. If individual 7 can send 
directly to 7 we say there is a link from i to j; this is an outgoing link for 1 
and an incoming link for 7. The structure of the communication net can be 
given by a matrix S with elements s,; where 
hie = Le, 

and fori 7 , 
Si; = 1 if there is a link fromi to 7, (2) 
$i; = 0 if there is no link from to 7. 


The units in the principal diagonal describe the fact that each individual 
remembers all the information he has received. 

We assume that S is connected, which means that for every ordered 
pair i, 7 there is a path from 7 to 7, that is, a set of non-zero elements 
Sii:SizigSizgig « « - Sqj- This is equivalent to stating that some power of S has 
no zero elements. The interpretation of connectedness is that completion 
of the task is possible, since there is a chain of individuals who can pass a 
message from any 7 to any 7. 

At each sending time, each individual chooses one recipient for his mes- 
sage, that is, one of his outgoing links. These choices give a sending pat- 
tern which can be specified by one of the set of sending matrices, S®, 

= 1,2,...m,. The elements of S® are s‘) where 


s@=1, 


and for each 7 

s() = 1 for one of the j’s for which s;3;=1, (7 4%), } (3) 
and 

s() = 0 for all other 7 #7. 


The assumption of random choice of outgoing links with equal prob- 
ability for each means the probability that 


5) = 1 (for i # 7) is ——*—— 


> s,,-1 
7 


since the denominator here is the number of 7’s outgoing links. The num- 
ber of possible sending matrices S® is 


m= T1( > su-1). - (4) 


Sij 


ities ee 
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Since the choices of the individuals are independent, the probability of 
any one of these is 1/,. 

The distribution of information among the individuals of the group at 
any time is specified by one of the information state matrices, C, with 
elements c‘%) where 


t7 


éeeil., if 7’ s initial information is known to 7 , 


(5) 


E{e) = Oise otherwise . 


At the start, the information state is given by C, the identity matrix; 
and at the completion time, the information state is given by the matrix 
whose elements are all unity—we call this C?. It is possible to write 
down 2"‘-) information state matrices, corresponding to the number of 
ways of placing 0’s and 1’s in the m(m — 1) off-diagonal positions. How- 
ever, not all of these can actually be attained because of the restrictions 
imposed by the S®. 

With this notation the information state which results if we have the 
state C® and send messages according to the matrix S” is the matrix 
product CS). In forming this product, Boolian arithmetic is used,* 
which means that 1 + 1 = 1, the other rules of addition and multiplica- 
tion being as usual. It is easy to see that the 7, jth element of Cs, 


which is 
de Ps 
k 


will be 1 if, and only if, there is at least one value of k for which ci? = 1 
and sf) = 1, but this means that ’s initial information is known to some 
individual, k, and & sends to j. 

The set of all the information state matrices, C“, is obtained by mul- 
tiplying C® by each of the S“; the resulting information state matrices 
are each multiplied by each S®, and this process is continued until no 
new information state matrices appear. 

We also need the quantity of information in the information state 
matrix C“. This is given by 

I= >), 


i,7=1 


* Boolian arithmetic is used only for these matrix products in this paper. 
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and simply counts the total number of pieces of information known to all 
the members of the group. 

4. Markov chain formulation. To treat the problem of calculating the 
RD we consider the C as the states of a finite Markov chain with sta- 
tionary transition probabilities (Feller, 1950). We denote the transition 
probabilities by @,g, this being the probability that when the information 
state is C™ it becomes C) after one sending time. These transition prob- 
abilities are easily calculated as follows: form the products C@S®, 


CMSQ,..., COS), if mg of these products are C®, then 
ng 
dag = —. 
oot a 


The matrix of transition probabilities (a.g) can be made to assume a 
particularly simple form, in fact, it can be written as a triangular matrix, 
i.e., all elements below the principal diagonal are zeros. This can be seen 
by considering J,. From the fact that s% = 1 for all 7 and vit follows 
that if C® = C™S©), so that 


cf) = p> c(2) so), 


we have c = 1 whenever c‘) = 1, so that Js = J.. Also a transition 
with Zg = J. can occur only if C®’ = C®. It is, of course, possible to 
have Ig = J, with C“) # C® but transitions between two such states 
can never occur. The only possible transitions are to the same state or 
one with a larger value of J,, so that by assigning the indices a in order of 
increasing J, (with any assignment for the C‘ having equal J,) the transi- 
tion probability matrix becomes triangular. 

It is now easy to characterize the states (Feller, Joc. cit.). The fact that 
structure is connected means that for every a except a = N there is a B 
with 6 >a and a.g > 0; for a = N, ayy = 1. Thus every state is tran- 
sient, except the final state (task completed) which is the single absorbing 
state. 

Let a{} be the probability of a transition from C®) to C® in ¢ steps, so 
the ai} are the elements of the ‘th power of the matrix, A = (a,,). Then 
a\? is the probability of completion of the task at time ¢ or earlier. The 
cumulative RD is then given by a‘, a2)... , and these values can be 


pas in succession by multiplying the row vectors (a, a, ..., 
ayy) by A, ice., 


t-F1 sy = 
Cagt) , GUF) yw SOLD) = (eG), alt), os gO) ay 
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However, we can do more than this. Explicit expressions for the gen- 
erating function of the RD and for the expected value of the completion 
time can be written. Let p,(#),¢= 1,2,...,anda=0,1,...,N—1 
be the probability that starting from C“? the final state C™ is reached in 
t steps and not fewer. Then 


Bol) Seno a= 0, 15 ee (7) 
N=1 

pa(t+1) = >> aspe(t), for t>0. (8) 
B=a 


Introducing the generating functions for the 7,(¢), 


Gilu= > Pali) ut, 
t=1 


we obtain using (7) and (8) 
1 foo) 
—Ga(u) = dan + 1 Patt 1) ut 


t=1 
aN! 


= den + >) Dy dashes (t) u! 


t=1 B=a 


= day + > >: be (t) ut 
B=a it 


so that 


“Ge (u) — 2 dapGp (U) = dan- (9) 


The set of equations in (9) can be solved in succession starting with the 


last ‘ 7 
(<- oy=1, 2) Gy (#) = Oy» 


(;- ayo, eo) Gy — “yee walwt (u) = Gy_5 (10) 


owes OS Ee ee ees eee MMe Li ut RC 
Ce 


1 9.) 6, (u) = 0yG, (H) ==, wiGiy = Boy 


GE t 
U 


We are mainly interested in G,(u), whose coefficients give the RD. 


G, (u) 
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From (9) or (10), we see that Gp(w) can be written in either of the two 
explicit forms: 


N-1 
aon oa dan 
Go(u) = + Da 
— —— (90 a=1 b= 200) iG a ace) 
U U (11) 
N—1 
eae Qoa Gap apn Ga, a+1 
aE 1 1 aaa | rete 
0<e<b ( = Goat Sa as8 on 6 ea 
U U U 
bd D9 Gos Bua ig aes Gon 
: ioe a 
yesh a5 a3 iY Ns 1N 
: a,, — a a 
1 0 ee a G3 24 2, N—1 2N 
pers 1 0 0 a ah a, — @ rr ale 
| (<- a 1 eh a a 3, N—1 3N 
U aa 
a=0 
1 
0 0 0 = @n-1, N-1 y @ y-1,N 


If E, is the expected value of the completion time starting from C“, so 
that E, = Gi(1), then from (9) using G,(1) = 1 and 


Sowa, 
B=a 
we obtain the equations for the £, 
peat | 
E.— >) apHp=1, a=0,1,...,N—1. (13) 
B=a 


Then E, = E(T), the expected completion time starting from the initial 
state C, is given by the expressions in (11) and (12) when w and all the 
a,y are replaced by 1. 

Equations for the second and higher moments can be obtained by tak- 
ing higher derivates of G,(a) using (9). 

5. Calculations. In Table I are given the generating functions of the 
RD and the average completion time for all possible connected nets for 
nm = 3 and a few for » = 4. 

The values of G)(w) and E(T) in the table were calculated using the 
formulas of Section 4, but in certain cases the work of calculation was 


(1 


TABLE I 


Structure Go(u) E(T) 
=I 
) uz 2 
FX 
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shortened because it was possible to represent the information state more 
simply than by using the matrices C“). This simplification is discussed 
below. 

We have also calculated, using (6), the first six values of p(t), for the 
linear chain with » = 5, i.e., the structure in Figure 3. These are com- 
pared with the values reported in CLM which were computed by a Monte 
Carlo technique on the Whirlwind I computer using 3000 trials. 


TABLE II 

t as 6 a 8 9 10 
polé)[Eq. (6)]........ 070. | 153 4) 0 Sige 2127 Be 
pod) (CLM). ...0... ‘078>| 140 | 197 | aso} 2119 Datos 


The calculations by these methods get very lengthy as 7 increases but 
the procedure is straightforward and could easily be handled by a high- 
speed computing machine, at least for moderate 1. 

For some special structures the calculation can be simplified because, 
for the present purpose of finding the completion time, the information 
state can be specified by less than the complete set of C“ as defined 
above. For example, in the star (Fig. 1) it can be seen that we need only 


ae 


ax 


FicureE 1 


specify the number a of individuals who have all the information. The 
distribution of completion times for this net is easily computed directly, 
but it can also be derived as a simple case of the Markov chain formula- 
tion. The result is 


n—2 


Go (u) past T (14) 


ae 


n—1 
E(T) =1+ (m1) S*=. (15) 
a=1 


~~ ew 


oa Te a}, NS 
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This can be easily generalized to the type of structure shown in Figure 2 
where there are the same number of points on each loop. If , is the num- 
ber of points on each loop not counting the central point and mu, the 
number of loops, then 

E SE 
(T) =2n,—- Vp Do 

Another example is the linear chain, Figure 3. The information state 
here can be specified by two numbers (a, 8) where a is the number of 
individuals who have received 1’s information and 6 is the number who 
have n’s information. This specification is sufficient because we are only 


f\ 
ea 
LS 


o<— 0 


FIGURE 2 


o<—>0<—>0<—>0°::+ 0<—>0 
| 2 3 nr 


FIGURE 3 


interested in completion time, so that if 2 sends to 3 before receiving 1’s 
information he must send again after receiving from 1. If, as mentioned 
below, we are interested in the time to reach an information state less than 
the completion state, this simplified specification may not be sufficient. 

A further simplification of the linear chain is possible because reversing 
the numbering of individuals gives an equivalent state, that is, (a, 6) can 
be identified with (8, a). In general, symmetries in the structure can be 
used to reduce the number of information states which need to be con- 
sidered. Thus in the structure where each individual has a link to each of 
the others, two information state matrices may be treated as identical if 
one results from the other by applying the same permutation to its rows 


and columns. 
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6. Generalizations and further problems. The assumption that each indi- 
vidual chooses from his outgoing links with equal probability can be re- 
placed by less restrictive assumptions without any essential change in the 
method. Thus instead of equal probability we could assign any fixed prob- 
abilities to each individual’s choices of links, or we could assign arbitrary 
probabilities to the sending matrices S, or the probability of any s™ or 
S®could even depend on the preceding information state C‘“. In any of 
these cases the transition probabilities a., are determined by C“™, and 
this is all that is required in order to have a simple Markov chain. 

However, the introduction of a rule such as: 7 does not repeat a message 
to j if it contains nothing new—does make an essential change. In this 
case d,g is not determined by a knowledge of C“ but depends also on how 
C“) was reached. Here we do not have a Markov process in the informa- 
tion states, C™. 

In some problems of communication in networks, the question is not the 
time to reach the complete state but a specified information state, say C*. 
This can be treated by the method of Section 4, if we identify with C* all 
the C“) which include C*, i.e., for which c) = 1 if cj; = 1. To insure the 
proper ordering, we must count the information content of C* as being 
n?, In this case, some of the simplifications of Section 5 may not be avail- 
able. The case where we do not start with C® being the identity matrix 
but with some other information state, C“ is also included above; in fact, 
we have the formulas for G,(u) and E£,(T). 

It was hoped that this method would give information about the struc- 
tures which minimize (7), but this problem will require further study. 
The values for small 2 which we have computed are probably not a very 
good indication of what happens as increases because of the rapidity 
with which the number of possible structures increases. However, an ex- 
tensive list of examples worked on a high-speed computing machine 
might indicate where to look for the best structures. 


APPENDIX 
Minimum completion time. The problem here is how to arrange the net 
and the sending schedule so as to achieve the shortest possible completion, 
T, subject only to the limitation that each individual sends not more than 
a fixed number, r, of messages at every sending time. We show that 


7 = smallest integer containing log,+:7 , (16) 
which is the same as stating that 7 is the integer satisfying 
log-4im S r< 1 +log yin 

nS (r+1)7< (r+1)n. j 


or (17) 


———— 1. +. 


eee. 
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It is easy to see that 7 can not be smaller than the value given by (16) as 
follows. Let J(¢) be the total information possessed by the group after 
sending time t. We have J(0) = x, and, if T is the completion time, [(T) = 
n?. Now after any sending time J(¢) can be at most r + 1 times its value at 
i eens PG) rr 1) (18) 
This is because the maximum amount of new information which each indi- 
vidual can give to the other members of the group by his 7 messages at 
time ¢ is r times the information he had at time ¢ — 1. Since this is true for 
every member of the group, J(#) < I(¢ — 1) + rl(¢ — 1), which is (18). 
The equality in (18) holds only if every message consists entirely of in- 


2 101 2 1 
id ol o————> 0 | fe) (@) 
mS wr 
2 
| ! | | | 00) 0 nee = SOIC 
(ome) 00 0o<——————> 0 10 
2 \ \ \ I 
000,0 <> 010 


vA SN 


100 2 ite) 


Ficure 4. Communication nets giving minimum completion time for n = 2, 4, 8, andr = 1 


formation which is new to the recipient; if there is any repetition of infor- 
mation already known then the inequality holds. From (18), 


I(t) S (r+1) #1 (0) =n(r+1)! (19) 


and using 7(T) = n? in (19), 

T 2 log,+1" 
and 

7 = log,+i" . (20) 

Equation (16) says that the equality in (20) is attained if 1 is a power 

of r + 1. To see this consider first the case r = 1 and let » = 2™. In Fig- 
ure 4 we have graphed the communication nets which give the minimum 
completion time for m = 2, 4, and 8. The numbers on the lines connecting 
the points show the paths over which messages are sent at each time. It is 
clear that the indicated pairing-off arrangement can be continued to any 
power of 2. Formally what we do is give each individual a number from 
0 to m — 1 and write the number in binary notation, i.e., 7 is written as 


e et (21) 


im ae m—1 i2 “a1 
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where each ¢ is either 0 or 1, and (21) means 


t= Cat2ent2est...+2™ 1eim - (22) 


The = 2” individuals require all the binary numbers that can be 
formed using m binary digits. The instructions for sending messages so as 
to attain the minimum solution time then are that at sending time ¢ each 
individual sends to that individual whose binary number is the same as 
his, except in the tth digit, that is, 7 sends to that 7 for which 


Cik = Crk for k=l, 24 et Lee eee 
(23) 


Cit a Cjt . 


It can now be seen that in at least m steps, the initial information of 
any 7 will be transmitted to any j. Suppose 7’s binary number differs from 
7’s in the mth, meth, ... , mth digits 


Cim, ~ Cjmy ’ Cim, ~ Cjmys +++ 9 Cim, ~ Cimz > 


and i and j have the same digits in the other positions. Then at time m, 7 
sends to 7, having the same binary number as 7 except that @;,m, = €jm,} 
at time m2, 7, sends to 72 with the same binary number as 7; except @i.m, = 
€jm,; continuing in this way 7’s information reaches 7 at time m,. This 
proves the equality in (17) for r = 1 and » a power of 2. 

Now if 7 = 2™ — 1, we use binary numbers as above except one num- 
ber is omitted, say 0. The sending instructions are the same, except that 
when the recipient is 0 according to (23) the instructions are modified so 
that now the recipient is the individual to whom 0 would have sent at the 
next time by the rule of (23). Thus at ¢ = 1, 1 sends to 10; at ¢ = 2, 10 
sends to 100; at ¢ = 3, 100 sends to 1000, etc. Thus all messages passing 
through 0 arrive at their proper destination one time unit earlier than 
before (when » = 2”), and since the recipient is the same individual who 
by (23) would send to 0 at the next time, all the information previously 
passing through 0 is passed to the entire group without loss. If» = 2™ — 2 
we omit the next binary number which never sends directly to 0, ice., 11 
in binary notation, and use the same modified rule for sending messages 
which by (23) would go to 11. This can be continued for all 2 between 
2™— 1 and 2™~!, the numbers successively omitted being 0, 11, 101, 110, 
1001, 1010, 1100, 1111, etc.; each of these is the smallest number which 


differs by two digits from all the preceding ones. This proves (16) for 
= 1, 
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For r > 1, the procedure is similar using (r + 1)-ary numbers. If 
n= (r+ 1), 


1 is written as eé. 


‘im 


€ ae. ey (24) 


“$,m—1°° v2 41 
where each e¢ is one of the integers 0, 1, 2, ... , y and (24) is interpreted as 
Oe ee (7 1) en et (Fe 1) ol es 2) 


The sending instructions are the same as (23) but now there are r num- 
bers satisfying (23). The reasoning showing that a message will pass be- 
tween every pair in at least m steps is similar to that for r = 1. 

For n = (r + 1)™ — 1 we omit the number 0, and modify the sending 
so that any message for 0 according to (23) goes at time 1 to any one of the 
number 10, 20, ... , 70 [in (x + 1)-ary notation]. It does not matter which 
of these is selected because each sends to all the others at time 2. Similarly 
at time 2, messages from 10, 20,. . . , r0 directed to 0 by (23) are now sent 
instead to any one of the numbers 100, 200, . . . , 700; these messages con- 
tain the information directed to 0 at time 1, and are included in the mes- 
sages sent to each of the numbers 100, 200, ..., r00 at time 3 by (23). 
Thus no information is lost. 

The numbers to be omitted for successively smaller values of 7 are 0, 1, 
2,...,7 — 1. Next we omit in succession 10, 20,..., (r — 1) O but re- 
tain r0. The next numbers to be omitted are 11, 21,..., (7 — 1) 1 but rl 
is retained. The scheme is that for each number omitted we retain at least 
one of its recipients at each ¢ according to (23). In this way any message 
directed to a missing number by (23) has a destination when the next digit 
is looked at in accordance with the modified instructions. All messages get 
passed on with no loss of information. It is also easy to see that the rules 
about which numbers must be retained allow us to reach any 7 between 
(r+ 1)™—1 and (r+ 1)". 

When 1 is not a power of r + 1 the sending pattern is not unique, the 
rules above being just one possible scheme. 
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The effects of structural features on various properties of enzyme systems are studied. 
Some of the effects are: in a homogeneous reaction, enzyme compartmentalization decreases 
the rate; in a heterogeneous reaction, compartmentalization increases the rate. The steady- 
state concentration of intermediates is larger in a non-uniform than in uniform systems. Perio- 
dicities do not generally occur in the common kinetic systems; they do occur in autocatalytic 
systems, but compartmentalization reduces their probability of occurrence. The conditions 
for overshoot are different for uniform and non-uniform systems. Multiple stable steady states 
are not a common occurrence among biologically typical reactions; they do occur in combined 
autocatalytic and surface systems (a mechanism for the gene position effect is suggested by 
this property). The local pH is affected by the enzyme aggregation as well as by the geometry 
of the enzyme structure. A 2-step system can give rise to the characteristic rate vs. pH curve, 
where the optimum is not necessarily at isoelectric point. The expression for the osmotic pres- 
sure inside a spherical particle is deduced. The pressure is shown to be dependent on the 
radius. The rate inside a cell particle is shown to be determined by the shape of the particle. 


Biologists have recognized for some time that a living system (cell, 
organ, organism) is not simply a homogeneous volume in which chemical 
reactions proceed. It represents, rather, a structure of definite spatial or- 
ganization, which organization is of decisive importance in determining 
the behavior of the cell. This fact has received additional recognition in 
view of recent developments in cytochemistry, where it has been shown 
that mitochondria are the seat of whole enzyme systems in a large variety 
of organisms, ranging from bacteria (Mudd, 1953) to insects (Levenbook, 
1953), plants (Bonner and Millerd, 1953), and mammals (Hogeboom and 
Schneider, 1950; Green, 1952). In particular, the mitochondria contain the 
enzymes of the all important tricarboxylic acid cycle (Kennedy and 
Lehninger, 1949) as well as of oxidative phosphorylation (Potter, Lyle, 
and Schneider, 1951), and the activity of these enzymes depends decisively 
on the spatial integrity of the mitochondria (Green, loc. cit.). 

Additional evidence concerning the existence of structural organization 
in cellular processes can be found in various areas of biology. We can men- 
tion here, among others, the localization of ribose nucleic acid in micro- 
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somes (Holter, 1952) which may imply the localization of protein synthe- 
sis. The localization in cellular membranes of enzymes participating in the 
transport of metabolites is also noteworthy (Lindberg, 1950; Rothstein, 
Meier, and Hurwitz, 1951). 

Another aspect worthy of mention is the discovery (Harman, 1950) of 
a definite correlation between the changes in the shape of mitochondria 
and alterations in the general physiological state of the cell. Here the old 
problem of the relation of structure to function has reappeared on the 
microscopic level, linked most intimately with some of the most impor- 
tant metabolic processes. 

It is interesting to consider the possible effects such structures may have 
on the behavior of the cell. In certain specific cases, as in the preservation 
of bioelectric potential or in active transport (Bierman, 1953), the func- 
tion of the structural elements can be fairly easily postulated. But in deal- 
ing with more general metabolic systems other, more general effects need 
to be considered. In a recent stimulating book, D. Danielli (1950) dis- 
cussed some of the variables which are of significance in controlling en- 
zyme systems. He listed the following factors: pH, sulfhydryl (SH) con- 
centration, concentrations of inhibitors and activators, and those processes 
controlling access of substates to enzymes. 

In this paper the following problems will be considered: 


1. The effect of enzyme localization upon the steady-state rate of a 
reaction sequence. 

2. The effect of enzyme localization on concentration of reaction inter- 
mediates. 

3. The influence of enzyme localization on competing substrates. 

4. The effect of enzyme localization on the manner by which the steady 
state is reached, e.g., oscillations and overshoot. 

5. The effect of structures on the possibility of multiple steady states. 

6. The effect of enzyme localization of the local pH. 

7. The effect of enzyme localization on osmotic pressure. 

8. Some relation between the shape of the cell particle and the steady- 
state rate. 


It is hoped that the consideration of these problems will lead to a more 
concrete evaluation of the role of structured enzyme systems. 
Before proceeding, we must specify what shall be considered a struc- 


tured enzyme system. If we consider cellular particles (mitochondria) two 
possibilities arise: 


. a. 


; 
7 
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1. The enzymes can be distributed uniformly inside the particle (most 
probably in soluble form). 

2. The enzymes can be fixed in specific regions or compartments of the 
particle. 


For our purposes the significant difference lies in the following. In the 
first case, the “uniform” system, each step of a reaction sequence is carried 
out in the same volume and hence diffusion is ignored. In the second case, 
the “non-uniform” system, different steps occur in different regions and 
hence products of one reaction must diffuse from the region of their pro- 
duction to the region of their consumption for the reaction to be able to 
proceed. On the basis of present experimental evidence (Palade, 1953) 
both types seem to exist. If the enzymes are uniformly distributed, this 
type of localization differs from the non-localized system (e.g., the en- 
zymes of anaerobic glycolysis) in that the reaction is compressed into a 
small volume and the intermediates are able to diffuse into and out of the 
reaction volume. The first model of a structured system is, therefore, a 
system in which all enzymes are uniformly distributed, but in which inter- 
mediates diffuse. Attention must be paid here to the problem of volume 
dependence of reaction rates. 

The second, “non-uniform,” model is that of a volume divided up into 
various compartments, i.e., regions in which enzymes are localized. As 
described before, for a reaction sequence to occur, substrates must diffuse 
from one such compartment to the next. This picture of a sequence of com- 
partments conforms to recent photographic evidence concerning the inter- 
nal structure of mitochondria (Palade, Joc. cit.). Electron microscope in- 
vestigation has shown that mitochondria are in all probability divided up 
lengthwise into segments or cristae, which could very well be loci of 
enzymes. 

These compartments may be so arranged that each enzyme is adjacent 
to the one following it in the reaction sequence, and, furthermore, are 
close enough that each product is simply “handed over” from one cata- 
lytic surface to the other. This is the picture biochemists frequently have 
in mind when they speak of enzyme localization (Michaelis, 1951). For 
these reactions, the advantage of enzyme localization is apparent. On the 
other hand, it requires a structure extremely well adapted to one particu- 
lar reaction sequence; that is, if the reaction proceeds in the order of inter- 
mediates, Xi, X2,.... Xn, the enzymes have to be located precisely in 
that order. 

We shall consider the possibility that the enzymes are not quite so pre- 
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cisely ordered, and that they are separated by some distance over which 
the substrates must diffuse. Since the distances involved are quite small, 
on the order of hundreds of Angstroms, a fairly random arrangement of 
compartments would be quite acceptable. 

We shall now proceed to examine the criteria listed above to determine 
the effect of enzyme localization on cellular metabolism and behavior. 

I. Steady-state rates. It is clear that one of the most significant properties 
to be examined is the steady-state rate of biological reactions. Such 
studies have been carried out for open and closed systems (Hearon, 1949a, 
b), but always on the assumption that the enzymes are uniformly dis- 
tributed throughout the reaction volume. We shall consider here the par- 
ticular effect of enzyme localization on the steady-state rate. 

Two comments need to be made concerning the model which we shall 
investigate. First of all, the volumes in which the reactions take place 
(mitochondria, etc.) occur in the cytoplasm, and it is this cytoplasm which 
makes up their external environment. We shall treat this environment as 
uniform and unchanging in time. Further, wherever experiments are sug- 
gested involving manipulation of the “external’’ concentrations of metabo- 
lites, we have in mind particles that have been separated from the cellular 
matrix by centrifugation and whose environment is now directly accessible 
to experimental control. Secondly, throughout this first section we shall 
treat reactions of the type 


k: 
Cit XsaBet+X;, F=1,2,...0. (1) 
me 


The C,’s, X;’s, X;1’s, and B,’s are all considered as substrate and product 
molecules. The enzymatic effect is included in the coefficients k; and k_;. 
This approximation is, of course, only valid as long as the enzymes are not 
saturated, i.e., as long as the reaction rate is linear with respect to the 
enzyme concentration. Hence, if so desired, we can also think of k; and 
k_; as v;H; and v_;H_; respectively, where v; and v_; are constants and 
E; and E_; the concentrations of the enzymes of the jth step. 

A. Uniform system. The first model, that of a uniformly distributed en- 
zyme system inside a cell particle, differs only from the commonly de- 
scribed open system reactions in that intermediates are permitted to dif- 
fuse. This case was treated very briefly by Hearon (loc. cit.), and we shall 
investigate it here more fully. 

Consider a reaction system as described in (1), where the concentrations 
C; and B;, as well as Xp and X,, are considered to be constant with respect 
to time, while the concentrations of the intermediates X; through X,,_; are 


4 
‘ 
. 
a 
4 
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functions of time. Let G(X;) represent the flow rate of the jth intermediate 
per unit volume out of the reaction chamber. By definition, the rates per 
unit volume of the jth and (j + 1)th step are given by 
0; = k,C;Xj-1 — k-;B;X; , 
. ; \ (2) 
U54-1 = R54i1Cy41X; = R—(54)Bj4iX j4i 

In the steady state dX;/dt must equal zero. Hence 

R3C5X 3-1 — R-;BjX5— heya y41 X53 + R-GavBj4iXj41—G(X;) =0. (3) 
It follows from (2) and (3) that 

Vj41 = 0; -G(X;), 

or, more generally, 


j+m—1 


Us4-@ UF a G(X;,). (4) 
k=7 


The meaning of (4) is clear. Given a system from which intermediates dif- 
fuse, different steps in the reaction sequence will have different steady-state 
rates. The relation between the steady-state rate of the (j + m)th and the 
jth step is given by (4). 

Our problem now is to investigate the effect of this diffusion of inter- 
mediates on the steady-state rate. We shall limit our investigation to the 
2-step system, because of the mathematical complexity involved in the 
solution of equations representing higher step systems. The system 
studied here is the same as that investigated by previous workers in the 
field (Hearon, loc. cit.) and can be represented as follows: 


Cit XoH— X14 Bi ) ; 
X1itQ,=—X2th. 


(S) 


Following our previous notation let w; = kiCi, w-i = k-Bi. As stated be- 
fore, X, is the only metabolite whose concentration is a function of time. 
In particular, we shall assume that it is kept at a fixed concentration X19 
outside of the reaction system, and diffuses out of the system at the rate 
(Xi — Xio) per unit volume. Consider a reaction chamber of volume 2V, 
surface area 2A, and concentration X\. If D is the diffusion coefficient of 
X,, the outflow per unit volume is given by 


2AD 


TVA, *1~ 4)» 


where Al, is an element of length of the order of magnitude of a particle 
diameter (Rashevsky, 1948). Hence \ = DA /VAh. We shall discuss this 
more fully later on. 
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We can now find the value of the steady-state rate of the first step, 
0, = WX» — w_1X1. Setting (3) equal to zero and remembering G(X1) = 
A(X = X10), we find 


_wXot AX +w-2.X2 (6) 
ee Wy+W2+A 
Inserting (6) into the expression for 2; gives us 
: _ Xo (A+ Wr) — w-\w_2X» — W-1AX10 (7) 
; A+. +W-1 ; 


If we designate by g the corresponding rate in the non-diffusing system 
(\ = 0), the ratio g/v, becomes 


aoe ieee wid [X10 (We + w_-1) — wW-2X2—W1X 0] (8) 
V1 (We + w_1) [wWiXo (A+ We) —w_yw_-2Xe —w—1AX10] © 


We shall also assume that in the non-structured system the 2-step chemi- 
cal reaction derived above proceeds in the forward direction. This requires 
that the change in free energy of that system be negative. This means that 
w—1W_2X2 


RT |n en 


is negative, or wyW2Xo > w4yw_eXe. It can be seen from (8) therefore that 
q/v is less than 1 for Xi) = 0, becomes 1 if 


W,Xo+w_2Xo 


Xio= 
Wo +w_ 


b] 


and approaches infinity at 


W Xo (A+ We) = Wwe aXe 
Vv iA 


X= 


For larger values of X10, g/v is always less than 1 (Appendix 1A). Hence 
variations in the external intermediate concentration result in wide vari- 
ations in the value of the above ratio. 

From (4) and (6) and remembering the definition of G(X), we can 
solve for v. The expression is 


_ WiW,Xo—w —jW»X_— AwWw_pX2e+ dw X10 | 


A+m+H4 a 
Hence q/v2 (q = % when \ = 0) is given by 
ie een Wer [wy Xo + w-2X2— Xro (we + w-y) ] 10 
Ve (w. ++ w—1) [wwe Xo + Aw2X io — W_-2Xo (W4+ A)” iy 


A number of observations can be made regarding (10). 
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1. If Xi is zero, and 


WW X 9 
W—» AG 2 


atin as (11) 


the ratio is greater than 1. Considering that 


W_\W_2 XG 2 


AF =RT In (12) 


WiWoXo ° 


the second condition can also be written \ < w_y(e~@7/8” — 1). 

If this second condition is not fulfilled, i.e., if the flow coefficient ) is 
large, this ratio becomes less than 1. 

2. The ratio becomes 1 when 


Wi Xo +w_.X2 


We + Ww_1 


Xp= 


and is less than 1 for larger Xio. It should be pointed out that g/v equals 
1 for the same value of Xi as g/2e. 

The effects of large Xi» and large \ are therefore the same: to increase 
V2 relative to gq. 

Systems such as described do, of course, occur in cellular particles, and 
it would seem possible, at least in principle, to measure the two rates, 2% 
and 2, independently (by measuring, for example, the rates of appearance 
of B; and Bz). The ratio of these two rates is given by 


as d [X10 (W-1 +2) — WwW Xo — w-2X2] (13) 
V1 WAXo+ XoWiwW2 — w_\w_»X2— w—1AX10 : 


This ratio can be plotted as a function of X10 (Fig. 1). Experimentally, 
four measurements should be possible: 

1. 2/2, when X49 equals zero, 

2. the value of Xi) when 72/7 equals 1, 

3. the value of Xi9 when 2 equals zero, and 

4. the asymptotic value of v%2/% as X10 becomes very large. 

If we assume that the second step is irreversible (w_. = 0) and desig- 
nate the experimental values of the four measurements above by Qi 
through Qu, the following ratios can be easily found 


W2 = Or 

r 1—Q,’ 

Ly ete 14 
= Um (14) 
a a ah 
; r ONO) : 
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Also, asa test of the validity of the model, the following relation must hold 


between the Q’s 
Qo = oe (15) 
B. The non-uniform system. We have defined the non-uniform system 


to be a reaction volume, with diffusible intermediates, which is distin- 
guished from the uniform system in that it consists of a number of com- 


Ficure 1 


partments, each of which is the locus of a particular enzyme. Further, as 
mentioned before, we shall assume these enzymes to be sufficiently far 
apart to make diffusion the principal agent of transfer from the site of one 
step to the site of the other. It was mentioned before that recent electron 
microscope observations of mitochondria (Palade, loc. cit.) show the exist- 
eee structures which could well serve as enzymatic loci (Schneider, 
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We will investigate the kinetics of such a two-compartment system in 
which the following reaction occurs: 


= — ky . i 
Ci+ Xo = Xu+8,, 
~ (16) 


X19 +f. én Lae By se Xoo . 


The first step can occur only in the first compartment because the enzyme 
catalyzing this step is localized there. The product of this reaction, Xu, 
diffuses out of this compartment, partly into the medium and partly into 
the second compartment. There, designated as X12, part of it reacts with 
the localized enzyme of the second step to form X22 and By and part also 
diffuses out into the medium. The effect of the enzymes catalyzing these 
reactions is expressed through the rate constants, k;, k_;. As in the uni- 
form system, we shall assume that Xi and Xy are time dependent vari- 
ables. For the purpose of this discussion Ci, C2, Bi, Bz, Xo, and Xz» are 
kept constant. In particular, we set ®; = kiC;, @-; = k-Bi. 

In line with previous notation, Gi(X1) represents the flow rate per unit 
volume of Xy into the medium, G2(Xy) the flow rate of Xu per unit vol- 
ume from compartment 1 into compartment 2, and G3(Xi2) the flow rate 
per unit volume of Xi. out of the second compartment into the medium. 
Hence the system is described by the following differential equations 


o¥tt = 9, X)—B4Xu— Gi (Xn) —G2,(Xu), 

(17) 
dXv _ 5 Xn—t Gs(X 
di = D_2Xo —MXwtG.(Xu) — 3 ( 12) 


We now set Gi(Xy) equal to (Xu — X10), Go(Xn) equal to Ao( Xu — X12), 
and G3(X12) equal to A3(X12 — X10). Since we assume that the volume of 
each compartment is exactly one-half of the total volume in the corre- 
sponding two-step uniform system, the flow coefficients in the two com- 
partments of area A and volume V can therefore be written as DA/VAls 
on the simplifying assumption that the two compartments are identical — 
with respect to permeability and diffusion properties. Furthermore, if we 
assume that the particles are cylindrical in shape and that the amount 
diffusing out through the ends is negligible as compared to that diffusing 
through the sides, Al, can be equated to Ah. Hence A; and 4; are each 
separately equal to of the uniform system. 

Setting (17) equal to zero, the steady-state rate per unit volume of the 
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first step, p: = m1 Xo — ®1Xu, can be easily found. The expression for 
X11 1s given by 


x Xo (W2+ Ne + Az) + AoW oXo2 + Xip (Q+ ArWe) (18) 
oe DW, (W_ it A+ Ag) +@_1(A2 +3) +2 : 
where Ve MA3 + MA a AoAz- Hence 


_ 0X0 [Q + We (Ar + Az) | — Xiow-1 (Ai®, +2) — Aww 1D_-2Xo9 (19) 
Q + We (Ar + Ao +W-1) +H-1(A2 + X3) ; 


Similarly, to find pr = WX 19 aa W_2X09, we first determine X12. It is 
given by 
Xp2= — 


DX ods + W_.Xo9 (Wi + Aart Ae) + X10 (Q+ A3W_-1) (20) 
Wy (W—1+ Ai + Ae) +W-1(A2t Az) +L : 


The expression for p2 follows immediately. It is 


__ WyW2X ode + W.X19 (Q+ A3W-1) — W_ Xo (Q+ W_-1dz +@_1)3) 


Za 
a Q + Be (Ar + Ao +W—1) + W1 (No + As) Ce 


These two results now enable us to compare the rates of uniform and non- 
uniform systems. To do this, we must of course compare rates for the 
whole reaction volume. This means that (19) and (21) must each be mul- 
tiplied by V, the volume of the reaction compartment, while the corre- 
sponding expressions for the uniform system, (7) and (9), must be multi- 
plied by 2V. Furthermore, some assumption must be made concerning the 
relative magnitudes of the concentrations of the constant metabolites in 
the two systems. 

As was indicated before, the only difference between the uniform and 
non-uniform systems is that in the latter the enzymes are concentrated in 
special regions inside the particle. Therefore there is no reason to assume 
that the concentrations Ci, C2, Bi, Be, Xo, and X» are affected by the 
enzyme localization. The only reasonable effect of this enzyme localization 
should be on the rate coefficients, k; and k_;, i.e., we should expect, as a first 
approximation, that concentrating the same amount of enzyme in half the vol- 
ume would double k; and k_;. Hence ®; = 20;, ®; = 2w_;. 

We can now evaluate the ratio 2Vv,/V 1, where 1 is the steady-state 
rate of the first step in the uniform system and /; the corresponding rate 
in the non-uniform system. We use the simplifying assumption that 
A: = 3 = AX. The following expression results: 


2% 4, W-1 (22+ rr) 
am i rare { 2 (ww Xo — w_yw_2 Xo) (22) 


+d (wi Xo — w-2X2) + AX io (We—w_y) }, 
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where 
Di = (A+B. + W_1) [201 X owe + 201 X or — 2w_yw_.X2 — 2w_yw_.Xo 


1 : (23) 
a 2w—1AX10 + 7s (2w wor X 0 + w,X or — 2W_—1AW>X 30 = W—1X 192) ] ‘ 


A number of observations can be made with respect to (22): 


1. If \. becomes infinite, the ratio of the steady-state rates approaches 
1. It should be pointed out that this is not a self-evident result because, 
even with infinitely rapid diffusion between the two chambers, the non- 
uniform system differs from the uniform one in that the compartmentali- 
zation of the two reaction steps is still preserved. 

2. The denominator D, is positive as long as fi, as given in (19), is posi- 
tive, i.e., as long as the external concentration Xj is small enough to per- 
mit the reaction to proceed forward. If this condition is met, and if all 
forward coefficients, w:X0, we, are greater than any of the reverse coeffi- 
cients, w_;, w2X», then the ratio of the steady-state rates is always greater 
than 1. 

3. If X, the flow coefficient to the medium, becomes very small, the 
ratio is always greater than 1. 


Since, in general, the forward coefficients do tend to be greater than the 
reverse ones, and since the flow coefficient is probably not too significant, 
we can conclude that the rate of the uniform system of the first step will 
tend to be greater than the corresponding non-uniform rate. 

A similar calculation can be carried through which compares the 
steady-state rates of the second step, i.e., 2Vv2/V pe. This ratio is given by 
2 V2 We (2wW4+ d) 


po uD. { 2 (wyw2X 0 — w_w_2X2) (24) 


+) (wyXo — w—2X2) + AX io (W2 — W-1) 5 
where 


Do = (A+ Wr +w_) [ Zeer X0+ 2 Aw. X19 — 2W-2X 2 (w_,+ d) 8) 


+o (2 hwow—1X 10 + W2X 10d? — 2w_1hw—2X 92 — A*wW-2X 22) | , 
2 


We can conclude from a consideration of (24) that 

1. the ratio approaches 1 for very large os, 

2. the ratio is always greater than 1 for small ), and 

3. if po, as given in (21), is positive, i.e., X10 bounded from below, then 
Dz is positive. If this condition is met, and if the forward coefficients are 
again greater than the reverse ones, the ratio is greater than 1. 

These conclusions, as can be seen, are very similar to the ones in the case 
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of the first step. We can therefore conclude that under certain assumptions 
concerning the order of magnitude of the reaction coefficients and flow 
coefficients, the steady-state rates of the uniform system will, in general, tend 
to be greater than of the non-uniform system. 

It may be rewarding to examine this conclusion a little more closely. 

There are two factors involved in the difference between uniform and 
non-uniform systems. One is that in the non-uniform system we posulate 
a finite flow coefficient; in the uniform system it becomes infinite. The 
second factor is the effect of compression on the chemical reaction. By de- 
creasing the volume in which a given particular reaction takes place, and 
hence by increasing the enzyme concentration, the rate per unit volume 
is accelerated. These two factors operate in opposing directions. It is the 
particular resolution of these two conflicting factors which determines the 
relative magnitudes of the steady-state rates in the uniform and non-uni- 
form systems. Since we have shown that in general the uniform rate will 
tend to be larger than the non-uniform one, we can conclude that the fac- 
tor of diffusion is more significant than the factor involving the compres- 
sion of the reacting enzymes into a smaller volume. 

That this explanation is correct can be seen immediately from the fol- 
lowing considerations. 

The steady-state rate of the general reaction system (1) can be written 
(Hearon, loc. cit.) 


% Tene BG [pr 
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ky =; y=, k_; =v_,E_; =0_; ; 
V represents here the volume of the particle, and N; and M; the total 
amounts of enzyme available. 

The total steady-state rate for the whole volume is then given by vV 
and is independent of the volume of the reaction system. Since the flow 
coefficient ) is also inversely related to the volume V, it can be seen upon 
inspection of (7) and (9), for example, that its introduction into the rate 
expression does not alter this conclusion. 

Since a decrease in the reaction volume does not increase the rate, the 
only effect of compartmentalization is to decrease the steady-state rate 
because of the finite flow coefficients. , 


( 
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At this point we note that the above results follow from a particular 
model of the reaction system, a model in which the essential acts, collisions 
between molecules, are assumed on the average to occur identically in 
each volume element. This is equivalent to the assumption that the reac- 
tion system is essentially homogeneous. This assumption is fairly adequate 
where the surfaces on which the reactions occur are relatively small as 
compared to the distances between molecules. On the other hand, when 
dealing with localized surface catalyzed systems it becomes less certain 
that this homogeneity actually holds. 

A different model, representing the non-homogeneous case, would be 
that of catalytic surfaces immersed in solution in which reactants and 
products diffuse to and from the active sites. Here the volume elements 
do not have the same properties. Diffusion occurs only in a certain region, 
free from chemical reaction; the chemical reaction occurs in a special area, 
on the catalytic surfaces. 

As the simplest example of this model, consider two surfaces, one at 
x = 0, the other at x = /. A substance A present in large quantities in the 
solution between 0 and / forms B on the surface at x = 0 at a rate k,A 
which we take to be constant. Substance B diffuses then to x = / where 
it forms irreversibly a substance C at a rate ks. We further assume that B 
can also re-form A at the first surface. 

In the steady state the problem becomes mathematically simply one of 
evaluating the one-dimensional Laplace equation with the boundary 
conditions 


-D(5-) Ge HOE 
dx 0 
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The solution is completely straightforward (Appendix IB). The rate of 
production of C is given by k,B(1), and hence equals 


Dkyk3A 
Dk3t+ Dkot koko” 


We can see that this rate is inversely dependent on /. As expected, the 


(27) 


(28) 


greater the distance between the catalytic surfaces, the slower the rate. 


Further, if the reaction is completely irreversible, i.e., 2 is zero, the 
steady-state rate is independent of the distance. 

Therefore in so far as the enzyme reactions inside of particles can be 
represented by such arrangements of surfaces, the distances between the 
enzymes become a significant factor in influencing the steady-state rate. 
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One way of looking at this effect is to say that the distance between en- 
zymes determines the “necessary” concentration of the substrate. In the 
above example, suppose we have two identical enzyme systems (two sur- 
faces each), differing only in the distance between the surfaces. To 
achieve the same rate, different concentrations of A are necessary in the 
two cases. If J, and & represent the respective distances, and A; and A» the 
respective necessary concentrations of A, we find, from (28), 


D (ha+ hs) + haksls nes 


Ata oe 


If the diffusion coefficient D of the substance B is large compared to 
(Roksl:)/(R3 + ke) this ratio becomes 1. But if diffusion is not so rapid, the 
necessary concentration Ai can decrease to Agl/Ie. 

We can also compare this surface reaction with the corresponding homo- 
geneous reaction. It can be shown very easily that if the distance / between 
the surfaces becomes small enough the surface rate will be larger than the 
homogeneous rate. 

In the previous example, let us assume that the number of active sites 
on a unit surface of that system corresponds to the enzyme concentration 
of a unit volume. Hence the steady-state rate, per unit volume of the 
homogeneous system, is given by 


kyk3A 
hat hy ca 
This is of course identical with (28) when D becomes infinite. Equation 
(30) must now be compared with the corresponding surface rate per unit 
volume, which in this case is simply #3B(/)/l. It can be seen that the 
smaller 7, the greater the surface rate. In particular, the two rates are 
equal when / equals 


— D(ko+hks) + V [D?(ko+ ks) 2+ 4 Dhoks (ko + ks) | 
For smaller /, the surface rate is larger than the volume rate. When D be- 
comes very large, (31) approaches 1, as expected. 

It can therefore be concluded that if we view the reactions from the 
standpoint of surface systems, compression of volume does result in in- 
crease of rate. And it can further be assumed that for localized systems 
such as mitochondria, in which the dimensions between enzymes are fairly 
small, the surface model is a better approximation than the volume model. 


See “es. 1. 
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II. Concentration of intermediates. Another criterion by which we can 
compare uniform and non-uniform systems refers to the total amount of 
intermediates stored in the reaction system in the steady state. This can 
be of importance biologically for the following reasons: 


1. At certain times intermediates may suddenly seep from the reaction 
system. Complete drainage makes it difficult to reestablish the desired 
steady-state rate. This seepage could be due, for example, to injury of the 
system, or to the necessity of adjusting to an unusual temporary environ- 
ment. 

2. Some of the metabolites reacting with the intermediates can also be 
used in other reactions. Hence competition exists for these metabolites, 
which will be greatly affected by the amounts of competing intermediates 
This follows from the law of mass action. 


Therefore for both of these reasons it seems useful for the metabolic sys- 
tem to have stores of intermediates available. 

Now consider the diffusing uniform system. The concentration of inter- 
mediate is given by (6). When \ becomes zero (6) becomes the correspond- 
ing concentration X, for the non-diffusing system. The ratio of the two 
expressions is given by 

es A [Xo (W2+w-1) — WX — W-2 Xo] 


== — ; oe 
oe a eres SECS ees sy (32) 


Clearly, if X19 is small, X;, the concentration in the closed system, is larger 
than Xj; for Xi greater than (wiXo + w_2X2)/(we + ws), X1 is larger. 
Hence, as in the case of the steady-state rate, the external concentration 
becomes a mechanism for adjustment of the intermediate concentration. 

We can similarly compare the concentration of intermediates in the uni- 
form and non-uniform systems. To make this comparison meaningful we 
have to compare total amounts of intermediate, which is given in the two- 
step case by the ratio 2X1/(Xu + Xw). This ratio can be computed very 
easily from (6), (18), and (20): 


Xut Xn 0; (33) 
X [2 (wyw,Xo — w_-yw_2X2) + X10 (w2—w-1) — A (W-2X2—-W1 Xo) ] , 
where 
Ds = (A+u2+w-1) Deka 4d weeX TED Xone . 

X [wrXo (2tdy-+ d) fw Xs (20-2) + AXio (r+ r+w-1)1 f. 
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A number of observations can be made with respect to (33): 


1. If \: becomes very large, the ratio approaches 1, as expected. 
2. If Nis very small and w, greater than wy, the ratio is less than 1. 
3. If w_1is greater than w» and ) is negligible, the ratio is greater than 1. 


Since on the basis of experimental evidence \ can be considered quite 
small, and the forward coefficients are generally considerably larger than 
the backward coefficients, we can conclude that the total amount of inter- 
mediate present in the non-uniform system will be larger than in the uni- 
form system. 

IIT. Access of substrates to enzymes. One of the factors listed by Danielli 
(Joc. cit.) as being significant in control of enzyme systems is the access of 
substrate to enzyme. There are two aspects to this: 

1. Itis advantageous to the cell that certain reactions proceed well. The 
cellular structures would then be so designed as to channel the substrate 
toward the desired enzymes. 

2. It is advantageous to the cell that certain reactions be suppressed. 
The cellular structures would then be so designed as to prevent access of 
substrate to the particular enzymes. 

It is clear on the basis of the previous discussion that a structure such 
as a mitochondrion can be very useful in focusing substrate flow toward 
the desirable enzymes and eliminating loss in seepage or undesirable side 
reactions. As a matter of fact, any system of membranes relatively im- 
permeable to the substrate and so arranged as to permit flow only toward 
the particular enzyme would accomplish this result. 

Another method whereby access of substrate to a particular enzyme is 
enhanced and access to an undesired side reaction prevented is spatial 
localization with respect to the substrate source. Consider the following 
simple situation. A substrate S is homogeneously dissolved in a medium. 
It can react with C; to form #1, or with C2 to form fs. The rate of the first 
step, mi, is kiC,S; that of the second step, qo, is keC2S. Hence 


(34) 


As seen from (34), in the case where spatial localization does not enter, 
the only factor determining the amount of substrate flowing to the desired 
goal, pi, and the undesired goal, po, is the ratio of the concentrations of 
Cy and ai 

Now consider the following localized system. As shown in Figure 2, the 
source of S is at A. It diffuses in from the source along the X-axis into a 
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compartment where only C; is present, or, if preferred, where the enzyme 
catalyzing the first reaction is localized. Hence, in this compartment, 


NY wl Os ee 21 ) 


where S, is the concentration of S in that region. Whatever remains of Si 
diffuses further toward B where the second enzyme is localized. Here p» 
is produced according to the reaction 


S2+C2— pro , 


where S> is the substrate concentration in the second region. It is clear that 


FIGURE 2 


such a spatial arrangement greatly enhances the accessibility of S to Ci 
and reduces it to Cs. The differential equations of this system are 


*- Ai (So — Si) — RiSiCi — Ae (Si—S2) , 
(35) 


——— de (Si a So) ar RoS2C2 . 


The ratio of the steady-state rates, @ = MCiSi, @ = keC2So, is found by 
solving for S; and S2 from the steady-state expressions corresponding to 


(35). It follows that ap ee 
oh 101 a2 36 
(1 + de ) oe 


To compare (36) with (34), assume that a certain ratio (q:/q2)* is desired. 
To make this feasible in the uniform system, a concentration 


_ HC: (2) Picts 
Gas ki \Q : 
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is required. Clearly, the larger C2, the greater C; must be. But in the lo- 
calized system, using (36), Ci is required to be 


(a sy" (“ie iG) (38) 


As expected, the smaller dz, the smaller the required Ci. Further, as Ce 
increases, C, does not become arbitrarily large, but reaches the asymptotic 


value of 
(2 
g2 


Enzyme localization, as described, is therefore a simple and effective 
mechanism for controlling reaction systems, enhancing certain reactions, 
and suppressing others. Such a mechanism could conceivably be involved 
in the well known transformation principle (Braun, 1953) or the alteration 
in bacterial synthesis upon viral infection (Epstein, 1953). In both cases, 
outside agents consisting largely of desoxy-ribose-nucleic acid (DNA) 
shift the synthetic pathways into new channels. It is suggested that DNA 
might be able to accomplish this, perhaps, by locating itself near the 
sources of energy and/or substrates being competed for, thus suppressing 
the normal enzyme pathways. 

We have mentioned so far two structural features involved in substrate 
accessibility, namely, permeability features (membranes, etc.), which pre- 
vent flow in undesired directions, and enzyme localization of the kind de- 
scribed in the model above. There is another way in which enzyme locali- 
zation could conceivably effect access of substrate to enzyme. It is well 
known that the addition of a colloid to a solution increases the solution 
viscosity. A particular expression for this viscosity, due to Einstein, is the 
well known expression 


7 = no (1+2.5¢), (39) 


where 7 is the coefficient of viscosity of the suspension, 7 that of the dis- 
persion medium, and ¢ the relative volume concentration of the suspen- 
sion (Gortner and Gortner, 1950). Since the diffusion coefficient is general- 
ly considered to be inversely related to the viscosity of the suspension in 
which diffusion takes place, it can therefore be concluded that enzyme 
localization could give rise to a considerable decrease in the rate of diffu- 
sion of substrate in the localization region. 

IV. Oscillations and overshoot. Although most of the work done so ate on 
cellular metabolic systems has dealt with the steady state, nevertheless a 
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good number of investigations have also been made of the time dependent 
behavior of these systems. Most attention has been paid here to the possi- 
bility of periodic behavior, and to the possibility of overshoot. 

A. Periodic behavior. Burton (1939), in his pioneering article, indicated 
interest in possible periodic behavior of open systems, but did not actually 
study it. A detailed investigation of the periodic behavior of autocatalytic 
systems was made by M. Moore (1949). She showed that those systems 
possess periodicities which are physically feasible. Hearon (1953) studied 
the possibility of periodicities in open and closed linear systems and came 
to the conclusions that, under certain general conditions, periodicities are 
not possible for them. 

We will attempt here to carry this investigation somewhat further, em- 
phasizing in particular the effect of localization on possible periodicity. 
Unfortunately no general, easily utilizable criterion is known to the author 
by which periodic behavior can be definitely established. 

The following procedure will be followed. In all cases dealing with vari- 
ous nonlinear systems, we shall assume that the system fluctuates only 
slightly from a steady state and hence we will be able to transform the 
original nonlinear differential equations into linear differential equations. 
These linear equations involve the amplitudes of fluctuations from the 
steady state. 

The solution of this set of linear equations will be of the form 2u, 
exp(— Ait) and the A, will be determined as the roots of the auxiliary equa- 
tion of the form 


N+ ay + a2... +a, =0. (40) 


The solution will be oscillatory if any of the 4,’s are complex; our problem 
therefore becomes one of determining the existence of complex roots of 
(40) without having to solve it explicitly. Unfortunately, the author is 
aware of only a sufficient condition for the existence of complex roots 
(Turnbull, 1944). If aj < 2a. in (40), complex roots exist. Of course, the 
condition may not be satisfied and yet complex roots will occur. This con- 
dition is particularly useful to us, because, given the determinant of co- 
efficients from which (40) is derived, we can easily determine the two co- 
efficients of A" and A*~, a, and ax. 
Consider a set of » differential equations 
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On the assumption that y, = Zu; exp (—A#), (41) becomes 
of S> yj (ass + 089) go Wh, Die AG OD 
ZAI 


where 5,; is Kronecker’s delta. When the determinant of the coefficients 
of (42) is expanded, the auxiliary equation (40) results. The coefficient a 
of \" is given simply by 


n 
> Qi. 
i=1 


The coefficient a of \”-? consists of two series. The first series is 


1 n n 
2 
2 s ii 055 — Bs 
‘=4 1 


and the second 


ix ; 


i.e., the sum of all products of all diagonal terms subtracted from the first 
series. It is clear that (Zaii)? > Laiiaj; — 2Za2;. Hence aj can only be less 
than 2a» if the second series of a2 is negative, thus canceling the minus sign 
in front. This can only be true if at least one a;; has the opposite sign of 
a;;. If this condition is not satisfied, we have clear proof that the sufficient 
condition is not satisfied. 

This criterion can now be applied to various systems. The following 
do not meet the above requirements: 

1. the linear uniform system, 

X:utC;=Bi+X:, 


as expected from Hearon’s investigation, 

2. the same system as in 1., except that it is made nonlifiear by treating 
the C,’s and B,’s as functions of time, 

3. an open enzyme system of the Michaelis-Menton type with the mod- 
ification that its second step be reversible, e.g., 


ky 
Ex Si= Xs ¢=1,2,...” 


XieSintEy, 


(Si is the ith substrate, X; the ith enzyme-substrate complex), and 
4. the non-uniform system, analogous to 1. and 2. above. 
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From the four cases above it can be assumed that periodicity is not a 
widespread property of simple kinetic systems. 

The following systems do satisfy, or can satisfy, the sufficient condition 
for periodicity: 

1. The closed reversible enzyme system. This differs from 3. in that the 
equation of conservation of mass is used. If K is the total amount of mate- 
rial present in the system, 


K= >) (Si+- Xi) + Sat. (43) 
X 
Setting a; = X; — X; and B; = S:; — Si, where X; and §; are steady- 
state concentrations, and using (43) to eliminate S,, the determinant of 
coefficients of the a,’s and B,’s now has one diagonal pair with opposite 
signs (Appendix IVA), namely, 


Qia+1) = — (ky — 0-1) (Ea — X1); \ 
Q(n+1)+1 = oit o—1S2 : 


Hence if #; > o_1, the barest possibility of oscillation exists. Clearly, the 
greater , the number of steps in the system, the smaller the possibility 
that this one term can outweigh the other negative terms. 

But it can be easily shown that oscillations are mot possible in the most 
favorable case, i.e., 7 = 1. The characteristic equation of this system be- 
comes (Appendix IVB) 


w—A{ Ry (E,+S;) +o,;+hk-1+ 0-1 ie Se) } a 
+£,{ o-1k1 (Ey + Si + Se) + o-1k-1+ o1k;} = () = 


Oscillations cannot occur if b? — 4ac > 0, where 5 is the coefficient of 2, 
a = 1, andc is the constant term in (44). By simple rearrangement of the 
coefficients of (44) it is easily shown (Appendix IVB) that 0? — 4ac is 
always positive. 

We can therefore conclude that closed reversible systems of the type 
described will not oscillate. 

2. It was shown by Moore (Joc. cit.) that autocatalytic systems show 
oscillatory behavior. She investigated m-step systems of the type described 
here as diffusible uniform systems. We shall therefore investigate similar 
non-uniform systems. 

Consider 7 linearly arranged compartments, as shown in Figure 3. The 
metabolites are designated A,;, where Ax is the concentration of the Ath 
metabolite in the /Jth compartment. In all cases A;; is produced in the jth 
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compartment from A;-1,; autocatalytically, ie., at the rate kA j;A jat,;- 
The metabolite 4,; diffuses out of the compartment into the medium at 
the rate o;;4;; and into the next (j + 1)th compartment at the rate 
r;(A;; — Aj,j41). The only exception is in the first compartment, where 
Ax is produced at the rate KAy. The differential equations describing this 
system are 


Agee 
CPE Dsl oe ors a _— 
aoe 5A; A; i+v RAs 414 ptt, 541 05,5414 5, 541 
Gta Pas Bs 
Pale ah fn) may (45) 
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The solution of (45) is straightforward and the determinant of coefficients 
easily ascertained (Appendix IVC). If Ai; represents the steady-state 


FIGURE 3 


value of A;;, we can solve for the two coefficients of \"-! and "2. They 
are 


2(n—1) 
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(46) 


2(n—1) 
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Hence 


2(n—1) n—1 


n=} 
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As can be seen from (46), the sufficiency condition for periodicity is not 
always satisfied. Since in the case of uniform systems studied by Moore 
periodicity always occurs, we can immediately conclude that compart- 
mentalization militates against periodicity. 
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It can also be seen from (46) (Appendix IVC) that the sufficiency con- 
dition is satisfied if 


A A*| 


pot pods gk ad 
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(47) 


Hence the flow coefficient 7; must be bounded from above. It should be 
pointed out here that this non-uniform autocatalytic system reduces to 
the uniform system by setting Aj,;,1 equal to A,;, which, using (45), 
eliminates all terms involving 7;. 

B. Overshoot. Burton (Joc. cit.) discussed in some detail the biological 
significance of overshoot phenomena. Some of the examples cited by him 
as evincing this property are: positive overpotential in nerve, oxygen con- 
sumption in bacteria and Arbacia eggs, and the so-called “vagus escape.” 
This same property of open systems was also investigated by Denbigh, 
Hicks, and Page (1948). Overshoot is defined as occurring when a system 
in moving toward a steady state from one side crosses the steady-state 
value and takes on values on the other side of the steady state. 

Mathematically, it means the following. Let the displacement of the 
system from the steady state be given by Zu; exp (—yit). Overshoot oc- 
curs, then, whenever this series changes sign. In general, we can see that 
the larger the number of terms in this series, the greater the possibility 
that such a reversal can occur. In particular, since a non-uniform system 
either always has 27 — 1 or 2n terms, the non-uniform system has a larger 
number of parameters which can be adjusted so as to make overshoot pos- 
sible for a corresponding uniform system of m terms. 

As the simplest possible example consider a 2-step system with uniform 
enzyme distribution. If we set x = 41 — Xi, where X;, is the steady-state 
value, the following differential equation results 


—“~= — (wi+wn2)x. 


The solution is x = ve~-+™", and overshoot is clearly not possible. In 
the equivalent non-uniform system, 
C+ Xoo XutB.; 
XptQo— XotBo, 
we set x1 = Xn — Xn, x2 = Xv — Xv, and from the resulting differen- 
tial equations find (Appendix IVD) 
X1= Mie M+ pee! } 


X2= Bien + Boe! , 


(48) 
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where 


_ 2A+meTWat ve a wa (49) 


1,2 


If we set x10 equal to x,(0) and x20 equal to x2(0), we can then solve for 
the coefficients in (48). In particular, for x: we find 


_ X10 (A+ w-1— ¥2) — AX20 
ca Vien V2 ‘ 


(50) 
2 AX20 — X10 (A+ 4-1-1) 


a AYA oe 


In order to have overshoot x; must be zero for some finite time ¢*. Ac- 
cording to (48) this requires that 


=e exp (y1— 2) ¢*. 
Me 


Hence if y1 > Y2, —(u1)/(u2) > 1. The latter requirement, considering 
(49) and (50), is equivalent to (Appendix IVD) 


ty —w—1— 2072 > [402+ (we —w-1) 7] 7. (51) 
1 

Equation (51) is feasible only if (x20)/(x10) < 0. Following the argument 
of Burton in his article mentioned above, we now examine the conditions 
for which (x20)/(x10) < 0 holds. In particular, we require that x19 and 
x20 represent differences in steady-state values of Xn and Xj. These 
steady-state values are given by (18) and (20), with the simplification of 
setting M => 3 ==,(() 

Inspection of these equations discloses that the inequality (x20)/(x10) < 
0 cannot be satisfied for any variations in Xo, Xo, or in any of the chemical 
reaction constants @;, @_;. It is, however, satisfied by a change in the flow 
coefficient \. In particular, if \ goes to \ + AX, x10 becomes ®AA(B_1H_Xe 
— ®%_Xo) and yx20 becomes —#_yAN(@1H»X_ — ®H.Xo). This is 
completely different from the situation found by Burton for a uniform 
system. There overshoot was only possible for a variation in the reaction 
constants, and impossible for a variation in the diffusion coefficients. We 
can therefore conclude that non-uniformity in the kinetic system alters 
the requirements and conditions of overshoot. 

We can now summarize the results of this last section: periodicities are 
not common characteristics in simple kinetic systems, at least as deter- 
mined by the sufficiency condition. Their occurrence is reduced in non- 
uniform systems as compared to uniform ones. In particular, the flow 
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coefficient \, unique to the non-uniform system, has a sign which militates 
against periodicities. 

Owing to the larger number of terms in the solution of the equations of 
non-uniform systems, more parameters can be adjusted to make over- 
shoot possible. In comparing a 2-step non-uniform system with a similar 
uniform one as described by Burton, we find that the conditions for over- 
shoot are altered as a result of the non-uniformity. 

V. Multiple steady states. Another property which is of biological inter- 
est is that of multiple steady states. We mean by this the ability of a kinetic 
system to possess more than one stable steady state and to move from one 
to the other if sufficiently displaced. Most steady-state systems discussed 
in the literature (Burton, Joc. cit.; Denbigh, Hicks, and Page, loc. cit.; 
Bertalanffy, 1950) show only one stable steady state, and for these systems 
the final state is independent of the initial starting point and only depend- 
ent on the fixed parameters of the system. This condition has been de- 
scribed as “equifinality,” as was pointed out in a previous note (Bierman, 
1954); this is not a general property of steady-state systems, since, clear- 
ly, systems possessing multiple steady states will move toward one rather 
than another of the possible steady states dependent on the initial starting 
point upon displacement. 

This property is not unique, but occurs frequently in nonlinear systems 
(Minorski, 1947; Rashevsky, Joc. cit.). It is also of biological interest be- 
cause shifts from one level to another do occur and are of decisive impor- 
tance, e.g., in genetics (mutations), adaptations, or embryology. Therefore 
it seems useful to inquire what kinds of systems can give rise to multiple 
steady states and how structure relates to this property, if at all. 

We have used the graphical method in this investigation, because it is 
the simplest. For example, given two variables with respect to time, x(#) 
and y(t), whose first-order differential equations are nonlinear functions 
#(t) = Py(x, y), yd) = Polx, 9), Pi, y) = 0, and P,(x, y) = 0 are then 
plotted in the x, y plane. The intersections of these two curves represent 
steady states. The stability of these steady states must then be deter- 
mined by going back to the original differential equations and noticing the 
effect of displacement from the steady state on the variables « and y. 

Clearly this method is really only applicable to two variables. For three 
variables the method can be extended to a three-dimensional model, but 
this becomes complicated already. We have therefore limited ourselves 


here to considerations of two variables. 
The following systems were investigated and were shown mot to possess 


multiple steady states: 
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1. a 2-step linear uniform system, 

2. a 2-step uniform system made nonlinear by letting B, become a 
variable of time, 

3. the corresponding nonlinear non-uniform system, 

4, a reversible 2-step enzyme system, 

5. a 2-step system where A forms B via the Langmuir adsorption iso- 
therm k,A/1 + oA, and is similarly re-formed from B at the rate koB/1 + 
o.B, and 

6. a 2-step autocatalytic system. This was treated by Denbigh, Hicks, 
and Page (Joc. cit.) and was shown to possess one stable and one unstable 
steady state. 


From the list above it is apparent that the simpler nonlinear kinetic sys- 
tems, both uniform and non-uniform, do not possess multiple steady 
states. It was possible, though, to find one 2-step system possessing two 
stable steady states; this system involves autocatalysis as well as surface 
reactions. It can be described as follows: 

Let x be produced via a Langmuir adsorption isotherm at the rate 
kyx/1 + ox. It then can either diffuse out of the system at the rate kyx 
or react with another metabolite y to form a complex xy. This complex 
goes over to y via a surface reaction described by koxy/1 + oxy. This is 
also an autocatalytic reaction by y. The metabolite y can also diffuse out 
of the system at the rate kay. 

This system, which possesses multiple stable steady states, is then de- 
scribed by the following nonlinear differential equations 


52 
dy Roxy (923 


dic, St eewy aalee 


Before proceeding with an analysis of (52), it is of interest to point out, 
first of all, that neither a pure, autocatalytic system nora purely surface 
reaction system alone gives rise to multiple steady states. A system that 
combines both those features, however, can do so. Secondly, both sur- 
face kinetics and autocatalysis are significant and characteristic features 
of biological reaction kinetics, and are moving more and more into the 
forefront of investigation in such areas as virology, genetics, and enzy- 
mology. Thirdly, surface reaction systems are generally indicative of 


enzyme localization and hence associated with what we called non-uniform 
systems. 
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The graphical analysis of this autocatalytic and surface system is given 
in Figures 4, 5, and 6. In each case the various intersections of the two 
curves represent the steady-state positions in the «, y plane. The directions 
of the arrows indicate the stability of these steady states. If the arrows 
lead back to the steady state, it is stable, i.e., any displacement gives rise 
to compensatory processes. If the arrows lead away from the steady state, 


| 


x 


FicurE 4. The case: ki < ks 


Ve 


O 
Ficure 5. The case: ks <i < hz + ka(o1)/(o2) 
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it is unstable. For details, see Appendix VA. It should also be mentioned 
that one steady state of this system is given by « = y = 0, which is not 
given as an intersection but shown as a point in the figures. 

Three possibilities must now be considered: 

1. Case: ki < ks. This system is represented in Figure 4 and described 
in Appendix VA. As can be seen from the graph, there is only one stable 


20 22 24 26 28 


Ficure 6. The case: k} > ks + ko(o1)/(o2). Evaluated for the values o; = o2 = 1, hj = 101 
ko = 20, kg = ky = 1. : , 


steady state, namely, at the origin. This, therefore, represents a transient 
system. 

yy Case: ks < ki < ks + ko(o1)/(o2). Here, as shown in Figure 5 and 
described in Appendix VA, exists one unstable root at the origin and one 
stable root for positive values of « and y. 

3. Case: ki > Rs + kx(01)/(o2). As can be seen from Figure 6, this sys- 
tem has four steady states, two stable and two unstable. Hence, given this 
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set of parameters, multiple steady states arise for this system. The system 


will occupy a particular steady state according to the particular values it 
has on displacement. 


We can use such a mechanism to suggest an explanation of the gene po- 
sition effect. Let « and y be the concentrations of two enzymes produced 
at a certain chromosomal locus F and let x have as a precursor a metabo- 
lite produced by a gene G at another locus. Hence the parameter f; in the 
action scheme above can be equated to w,G, where G is the precursor con- 
centration at F. This value is given by the following equation: 

dG w,Gx 


may a. (Gy — G) Sis rs 


(53) 
where Gp is its concentration at the locus of G. In the steady state 
AGo 


Wx 
SST eres § 


Also, \ is given by Q// where Q is a constant and / the distance between G 
and F. Hence if the distance / is increased because of a mutation involving a 
change in gene positions, G is correspondingly decreased. But, as pointed 
out, &; is proportional to G. The enzyme system may therefore start out 
in state 3 with its two possible stable steady states and, as a result of a 
shift in the respective distances between the two genes, the value of k; may 
be sufficiently decreased to make only state 1 or 2 possible. The system 
must then leave its previous steady state and drop to a lower one, includ- 
ing, possibly, to the zero steady state. 

It can be concluded therefore that multiple steady states are possible, 
but require somewhat unusual kinetic systems. In particular, a combina- 
tion of autocatalysis and surface reactions gives rise to two stable steady 
states if certain inequalities are satisfied. Such a model can be used to ex- 
plain gene position effect. The required property of surface reaction again 
indicates the importance of the role of structural features in biology. 

VI. pH and structure. The significance of pH in cellular systems is well 
appreciated and a great deal of experimental and theoretical work has 
been done bearing on various aspects of this question. In this paper we 
shall investigate the following two problems: 

1. It is well known (Danielli, Joc. cit.) that the amphoteric nature of 
proteins (enzymes) gives rise to a local pH near the surface of the protein 
which can differ markedly from the average value in the medium. We 
shall study the effect of enzyme aggregation on this local pH. 
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2. It is also well known that enzyme reactions are very sensitive to 
changes in pH. In particular, we shall examine the effect of enzyme locali- 
zation in a 2-step system on the optimum pH for the reaction. 


A. The local pH. Since enzymes are amphoteric in nature, they can 
possess a net charge. This net charge attracts ions from the medium to 
form a diffuse double layer as described by the Debye-Hiickel theory. In 
particular, H+ and OH~ ions will be attracted, thus giving rise to a local 
pH capable of differing considerably from the average. 

The enzyme aggregation causing this electrostatic field can, of course, 
differ in extent. It seems reasonable to assume that the size of the aggre- 
gation will have an effect on the magnitude of this local pH. 

We shall investigate two models. The first is represented by a sphere 
whose outer surface consists of a layer of charge-bearing enzymes. If we 
increase the radius of the sphere and keep the surface charge density con- 
stant, we are essentially describing a situation of increasing enzyme aggre- 
gation. Physically, perhaps such a sphere could correspond to a microsome. 

We now assume that the active groups on the enzyme surface are react- 
ing with the Jocal H+ and OH~ ions according to the following scheme 


Ht +RA-s RAH,  OH>+R’Ct=R'CO+45,0, 


where RA™~ represents a surface anionic group, R’C+ a surface cationic 
group. It therefore follows, remembering (H*+)(OH-) = K = 10-4, that 


Shae 

(H+ (54) 
(R’CT) = a3 (H*)., 

where a; and az are constants. Physically, this means that the surface 

charges of the enzymes are effected by the local pH of the medium at the 

same time that the local pH is determined by the surface charges. 

Let «» be the hydrogen ion concentration in the medium, i.e., r= ~, 
and let ¥(r) be the electrostatic potential due to the charged enzyme sur- 
face. It follows then from the Boltzmann distribution that the hydrogen 
lon concentration at r = a, i.e., at the sphere surface, is given by 


X= Xp exp [—" i | 


(RA>) = 


’ 


hence 


(OH-] -*x Ens i ] : 


Here e represents the charge of one electron, # the Boltzmann constant, and 
T the absolute temperature. Remembering now that the surface anions 
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and cations react only with the H+ and OH™- ions at r = a, we can 
rewrite (54) 
ets Bn ns ep (a) 
(RA~) = = OP fe. | 
(55) 

(RCry = AoXq EXP | - ely . 
Our problem now becomes that of finding y(r). This can be determined 
from the Debye-Hiickel theory. It is given by (Verwey and Overbeck, 


1948) 
vy (7) = (ear (56) 


where K is the well known Debye-Hiickel constant. The coefficient Yo is 
evaluated by using the relation 


am (¢ Bu , (57) 

dr € 
where ¢ is the dielectric coefficient. The surface charge density o is given by 
g=—e€|(R'C*) — (RA). (58) 


From (55), (56), (57), and (58), the local pH can be determined quite 
simply. On the basis of the assumption underlying classical Debye-Hiickel 
theory, namely, [ey(r)|/(AT) «1, the resulting transcendental equation 
can be linearized. Our final expression for the local hydrogen ion concen- 
tration [H+](a) is given by (Appendix VIA) 


dare? a (a2%5 — a) (59) 
H+] (a) = i 
[H*] (a) = %o exp er | e(aK-+1) m+n a (asx + an) | 


It is clear from (59) that if x» > v/[(a1)/(a2)], the local hydrogen ion con- 
centration decreases with increasing radius a, and has the asymptotic 


value 
4a e (a2x5 — a1) (60) 
+ = = ‘ 
[H*] (a) = %o exp aT [eKxo+ Be paca a1) | 


On the other hand, if #» < V [(a1)/(ae)], the local hydrogen ion concentra- 
tion increases with increasing radius a, and has the asymptotic value 


4 e° (a1 — a0) (61) 
+E = : 
Ee at), setee SSP er [eKxot Fy (asxh+ a) | 
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Hence if the pH of the medium is great, the local pH will decrease with 
increasing enzyme aggregation; if the pH of the medium is small, the local 
pH will increase with increasing enzyme aggregation. These conclusions 
are summarized in Figure 7. 

The second model which we shall investigate differs from the first in 
that the enzymes are assumed to be adsorbed on the inner surface of a 
spherical shell. The significant pH is therefore that on the inside of this 
sphere. Physically, this may correspond to a structure such as a mito- 
chondrion. 


(H+), 


O 


a 


Ficure 7, Hydrogen ion concentration vs. radius; enzyme on outer surface 


The potential inside of such a sphere is given by Debye-Hiickel theory 
(Appendix VIB) as 


V(r) = sinh ar: (62) 
and the constant Yo is determined from 


dy _ 4no 
CTS Piet Wie 


(63) 
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From (55), (62), and (63), Yo can be determined. Again using the approxi- 
mation (e/kT) ¥(a) <1, we finally determined (Appendix VIB) 


By 


y (a) = i om 
(coth Ka—z-) K+ 8B, 
where 
= = (axt0— 1) 
Xo 
Rees (65) 
eS : Qi 
Bz = aT (a:t0-+*) 
Since 
ee ss Mande inky (nie 
Ae Bs oe K+ Bo’ 


we can see from (64) that (a) is monotonically increasing for 6; > 0 and 
monotonically decreasing for 8; < 0. Since the hydrogen ion concentration 
is given by 


Xo exp| —Za¥ (a) | : 
the following conclusions can be drawn: 
1. If x is less than +/[(a1)/(a2)](B: > 0), the local (H*) starts off, for 
a = 0, at 
e By 
Xo exp| = ry 3 | : 
which is greater than Xo, and decreases to 


é By ee 
Xo exp| — oT Kes for a= Oo, 


2. If xo is greater than ~/[(a1)/(a2)](6i > 0), the local (H+) starts off, 
for a = 0, at 
won 
aed GES 
which is less than 2, and increases to 


Xo exp| eral ; 


These conclusions are summarized graphically in Figure 8. 
The behavior of the local pH inside of a sphere is therefore considerably 


different from its behavior outside of a sphere. 
B. pH and a 2-step enzyme system. It is well known that changes in 


pH greatly affect the rates of enzyme reactions. This is easily understand- 
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able, since a good part of enzyme-substrate interaction is probably medi- 
ated via electrostatic forces which are very much affected by the electro- 
lyte concentration in the medium. Also, hydrogen and hydroxyl ions prob- 
ably react with charged groups on the enzyme surface and hence compete 
with the substrate for reactive sites. 

It has been shown for a good many enzyme systems that an optimum 
pH exists, differing of course from enzyme to enzyme, and that the reac- 
tion rate drops away quite sharply on each side of this optimum. Various 


(H+), 


Qa 


FicureE 8. Hydrogen ion concentration vs. radius; enzyme on inner surface 


theories have been advanced to account for these characteristic curves; 
some of them indicate that the expected optimum would be at the iso- 
electric point for the protein. Actually a good many of the experimental 
curves for enzymes place the optimum pH not at the isoelectric point but 
as much as 1 or 2 pH units off. This would seem to indicate that the rela- 
tion between pH and enzyme action is of a more complicated nature. 
We shall attempt to derive the typical bell shaped rate vs. pH curve 
from a simple model which does mot require that the optimum pH be at 
the enzyme isoelectric point. To simplify the calculations we shall assume 


7 
> 
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here that the hydrogen and hydroxy] concentrations are uniform through- 
out the whole region, i.e., the Boltzmann distribution according to the 
electrostatic field is ignored. Our model is as follows: Assume two enzy- 
matic surfaces placed a distance / apart, so that the first step of a reaction, 
S; — pi, occurs on the first surface. We assume that S; is present in large 
enough quantities in the medium so that its concentration can be consid- 
ered constant. The substrate p; diffuses through the medium and reacts 
on the second surface to form #». We shall also assume that the production 
of p; occurs on negatively charged sites of surface 1, while the production 
of p. occurs on the positively charged sites of surface 2. Hence hydrogen 
ions will act competitively on surface 1, and hydroxyl ions will act com- 
petitively on surface 2. We shall also assume that 9; tends to re-form 5; 
on the first surface, but via a different uncharged enzyme. 
The reactions taking place on the first surface are, therefore, 


Fy — ee (ES) ’ 


(ES;) 3p: +E., 
E,+H*s3 (EH), 
ke 
Pi Le Sy . 
The steady-state rate of production of /; is then given by 
o2 (ES;) — kop: (0). 


The first part of this expression, 02 (ES;), follows directly from the theory 
of competitive inhibition (Wilson, 1949). It is given by 


x a0. Ey r 66 
o2 (ESi) = (1+ KiH*) (o2+ 0-1) + 0181’ Se 


where E,r is the total concentration of enzyme available on the first 
surface. If p,(/) represents the concentration of p; at the second surface, 
we have a similar reaction there, giving rise to pe with OH~ ions competing 
for the positively charged active sites. The rate of p: production is there- 
fore given by 


o304Eorpi(l) _ 67 
vip (1) = 1+ K,0H-) (04 o-s) +0352 te 


To simplify calculations we shall assume that p,(J) is much smaller than 


+ (1+ K.0H-) (os 0-2) - 
3 
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Hence (67) can be approximated by 


a o3 04H ri (1) 68) 
vp: 1) = TK, 0H-) (oi ons) ( 


The desired rate is given by (68). But to evaluate it we must find #;(/). 
This requires solving 


ap, _ 
dx? 0, 


with the boundary conditions 


— 41) 9, ES) —beps(0),  — FE) = ap. 
These equations in p; are identical with (27) and (28) in B, except that 
ki A = o2 (ES), k3= 1. 

If we again set H+ = x, OH- = a/x, the steady-state rate, 1p,(/), using 

(29), can finally be written 


_ = ee (69) 
1 aE Yat + Yat? 


where 


V1 = 0102030491; 7E27D , 
2 = DkzK204(o4+ 0-3) (02+ 0-1+ 015i) , 
¥3a= (02+ o-1+ oS;) [Dk (04+ 0-3) + o304E27 (D+1k:) | (70) 
+ Ky (02+ 0-1) 04 Dk2 Ke (04+ 0-3), 
v4 = Ki (02+ o-1) [Dko (o4+ 0-3) + (D+1hs) o304E 27) . 


As can be seen from (69), v3 is zero for « = 0 and is again zero forx =~. 
In particular, the optimum hydrogen ion concentration is given by 


re V(2) a 


It is clear from (70) that this optimum is determined not only by the pa- 
rameters involving the enzyme itself, but also by the separation between 
the two enzymes, and the diffusion coefficient D of the substrate p,. Hence 
the optimum is not necessarily at the isoelectric point but is determined, 
in a complicated manner, by the total system. 

VII. Osmotic pressure. An integral part of studies on electrolytic solu- 
tion has been the investigation of osmotic pressure due to such phenome- 
na. One of the best known approaches to its determination has been via 
the Donnan equilibrium. 
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We shall examine this problem from our particular standpoint of en- 
zyme aggregation. Consider a sphere of radius a, whose inner surface is 
covered with protein material bearing a net charge. As a result of the po- 
tential arising from the surface charge, electrolytes from the solution will 
be attracted to or repulsed from the region of protein concentration de- 
pending upon their sign and the sign of the net charge. If the result of this 
attraction and repulsion is an increase in the total number of molecules in 
the inner spherical region, as compared to the medium, osmotic pressure 
will develop and water will tend to flow into the sphere. If the total num- 
ber of particles in the enzyme region is smaller, water will tend to flow out 
into the medium. As a first approximation, we shall use the van’t Hoff 
equation for osmotic pressure: 


P= >> (nz — Ni) ’ 
i=1 


where 7; represents the average number of molecules of the th electrolyte 
per unit volume of sphere, and min its number per unit volume in the 
medium. 

Further, 2,(r), 0 <r <a, is given by mo exp [—ez.b(r)/kt], where 2; 
is the valence of the ith electrolyte. Hence the average number per unit 
volume is given by 


pet 92 nig e-WAMOOVET Ar. (72) 

a’ 0 
This intergral can be evaluated by expanding the exponential and discard- 
ing all but the first three terms. Remembering (62), the final expression 


becomes (Appendix VIT) 


—— VJ (Z22nio) (sinh 2Ka—2Ka), (73) 
where K = d+/(222ni0). We observe from (73) that if Yo is independent of 
the radius P becomes infinitely large as a becomes infinite. 

We now modify the picture, by assuming further that the electrolytes 
in the medium not only form a double layer, but also react chemically with 
the charged groups on the enzyme surface; hence affecting the charge 
density and, therefore, Yo also. Using (63) we can, therefore, evaluate po. 
As shown in Appendix VII it becomes 


ee RE Ne (74) 
vo= aK cosh Ka+sinh Ka(aé— 1) 
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where 


dnl Br Ar ] cee | Br Ar | 
on ES: : ekT LQ;n-: DK M4: J 


€ 
We can now express P as a function of a. It is given as 


es Gs 
A BTN 


= a (sinh 2Ka— 2 Ka) 75 
P a/ (Zitio) [aK cosh Ka+sinh Ka (aé,— 1) 12 (75) 


A few observations can be made about (75): When a equals zero, P is 
zero. When a becomes very large, P becomes zero also; hence P must have 
a maximum for some radius a. 

It is difficult to assess how significant the actual magnitude of this os- 
motic pressure could become. It is of interest to read in the experimental 
literature (Shelton, 1953) that mitochondria, as observed in cultured 
fibroblasts, ‘‘assume a bewildering variety of shapes. . . . Large ‘hollow’ 
blebs appearing at the ends or in the middle of a mitochondrion would 
swell and collapse, frequently pinching off a piece of mitochondrion at the 
point of rupture.” 

It is tempting to speculate that perhaps these morphological changes 
could be due to variations in osmotic pressure, which variation might per- 
haps be traced to changes in charge density of mitochondrial enzymes as a 
result of a progressing reaction or of alterations in the size of the particles. 

It also seems possible to ascribe the lysis of viral hosts to osmotic pres- 
sure resulting from the production, by the virus, of a large number of high- 
ly charged DNA molecules inside the cell. The lytic action of phage 
would, therefore, be due to the act of localizing electrolytes inside the 
host cell, hence giving rise to an increasing charge density in that volume. 

VIII. Morphology. As previously mentioned, a number of investigators 
have found some correlation between shape of mitochondria and varia- 
tions in metabolism. For this reason it is of interest to see whether any 
simple conclusions follow from our kinetic analysis concerning relation be- 
tween shape of these cell particles and steady-state rates. 

Unfortunately the limitation of the kinetic analysis to 2-step systems 
makes any general deduction rather dubious. Also, there are many other 
factors in addition to simple steady-state rates which play significant roles 
in this relation and which are not considered here. 

Considering rate only, then, we see that the rate in a certain region will 
be very much determined by the amount of diffusion possible, and this, in 
turn, will depend on the total surface area available for outward diffusion. 
The sphere having the smallest surface for a given volume will have mini- 
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mum diffusion. In the case of a 1-step open system with one diffusible inter- 
mediate, this becomes the only significant effect and therefore the rate of 
the sphere will be minimum, as compared, for example, to the rate in a 
cylinder of equal volume. 

Consider now a 2-step system with two diffusible intermediates of the 
form 

BCE. 

Here the situation is somewhat different. If the surface area is too great, 
all of the first intermediate C, diffuses out of the system and the second 
step never occurs. Hence there exists an optimum surface area for a given 
volume, as distinct from the previous case where the optimum was the 
maximum surface area possible. 

This can be seen more concretely by using Landahl’s (1953) approxima- 
tion method. The rate of appearance of C2 in a sphere is given by (Appen- 


dix VIII) 
16 Ky Kor DR® 


Wi= Tae D2+4 12 DR®(K-1+ Kat Ka) + K-1K—2R*” eh) 


where D is the diffusion coefficient and R the radius of the sphere. In the 
case of a cylinder of volume V = $7R* = 7hj20, the rate of appearance of 
C; is given by (Appendix VIII) 


6.2 6 2 
Wp = OT DV na (4h VY 17) 
‘iP 
where 
P= [124Dh5+ 6DV?2+ K_,Vh?] [24Dh5+6DV?+ K_,+ K,V*h]] 
— K_,K,hiv*. 


Clearly W2 is zero for both i) = 0 and very large iy. Hence an optimum 
value of hp exists, as expected. Further, there will exist at least two ranges 
of ho for which the rate of the sphere is larger than that in the cylinder. It 
is also possible that the rate in the sphere is greater than that in the cylin- 
der for any value of fp. It is clear, therefore, that the conclusions drawn 
from a 1-step system are not applicable to higher order systems. 

IX. Discussion and conclusions. It has been the purpose of this paper 
to examine the significance of structure of enzyme systems. The biological 
importance of this problem can hardly be overestimated. It is becoming 
more and more apparent as a result of biological research that life means 
organization and pattern—not only chemically, not only in terms of a 
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given set of reactions, but also in terms of spatial and geometric relation, 
linking together chemical reactions, bioelectric potentials, mechanical 
forces, and all the many processes characteristic of living things. 

It is becoming apparent that these structures must be studied not only 
statically, although this is of undeniable importance, but dynamically, in 
their relation to the ongoing activities of the cell. This sentiment is well 
known today and needs little discussion. But it is not too clear how this 
problem should best be attacked, and most of our knowledge concerning 
the relation of organization to function is still on a far too vague level. 

This paper was conceived as an attempt at a particular approach to- 
ward this basic problem. In all cases attention was focused on how to re- 
late a specific biological, chemical, and physical process or property to a 
given structure, or changes in structure. The purpose is frankly explora- 
tory; the author does not believe that biology is at the stage where pene- 
trating theoretical developments can occur apart from an integrated oper- 
ation involving experimental and theoretical work. 

On the other hand, such theoretical work can open up lines of thought 
and indicate relations which otherwise might be missed. It also offers op- 
portunity for collecting diverse data and arranging them in a semblance of 
order. Further, it develops models and their consequences which can be 
of significance to the field. 

This paper was devoted to a series of questions, all focusing on the prob- 
lem of the relation of structure to activity. The first question concerned 
itself with the effect of enzyme compartmentalization on the steady state 
of rate of reactions. It appears that for the particular kinetic system 
studied, the delaying effect of diffusion is dominant, and therefore that the 
steady-state rate is decreased as a result of this kind of compartmentali- 
zation. This result led quite naturally to an examination of the validity of 
the kinetic scheme used. It then appeared that the conclusion was perhaps 
of dubious validity because the basic model was biologically inadequate. 
It was suggested that perhaps the more fitting model for reactions inside 
of cells, in particular inside of particles, would be surfaces immersed in 
solution. A very simple examination of this model indicated that here 
compartmentalization does lead to an increase in the steady-state rate. 

Another significant problem dealt with was that of access of substrate 
to enzyme. It was shown that a simple localization of enzymes with respect 
to the enzyme source can give rise to important effects on the competitive 
situation, and assist the channelization of reactions one way or the other. 
The effect of enzyme localization was also studied with respect to the 
steady-state concentration of intermediate, and it was shown that in gen- 
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eral the non-uniform concentration will be larger than the uniform, hence 
also assisting competition. 

We also dealt with the question of possible periodicities. It was estab- 
lished that the simple kind of kinetic systems, both uniform and non- 
uniform, do not satisfy the sufficient criterion for periodicity. Further, in 
the case of autocatalytic reactions, diffusion which enters as a factor in the 
compartmentalized system seems to reduce the probability of periodicity. 

The phenomenon of overshoot was also examined and it was shown that 
structural features, if involved in reaction systems, provide additional 
parameters which can lead to overshoot. Furthermore, it was shown that 
the conditions for overshoot in terms of variations in the system parame- 
ters are different for uniform and non-uniform systems. 

Some space was devoted to multiple steady states, a subject of biologi- 
cal interest which has been insufficiently investigated. It was shown that 
a number of systems which commonly occur do not possess this property, 
but that a combined surface and autocatalytic system does. A possible 
mechanism whereby such multiple steady states could give rise to gene 
position effects was postulated. 

The problem of pH was then investigated, again with particular focus 
on the relation between pH and structure. It was shown that the local pH. 
in the neighborhood of amphoteric enzyme systems is affected by the ex- 
tent of enzyme aggregation. The nature of this effect depends on whether 
the enzymes are spread on the inside or outside surface of the particle. 
Specific relations between changes in aggregation and local pH were 


established. 


Another well known problem which was investigated with reference to 
structure was the relation of reaction rate to pH and the nature of the op- 
timum pH. In particular, a 2-step surface system was studied in which the 
reactions were assumed to occur via differently charged groups. It was 
shown that the reaction rate, when plotted against the pH, showed the 
characteristic curve, and that the optimum pH was not at the isoelectric 
point, but determined by the parameters of the whole system. 

The problem of osmotic pressure effects resulting from charged surfaces 
was also examined. The expression for the osmotic pressure inside of a 
sphere of uniform charge density was deduced; it was shown that the os- 
motic pressure has an optimum radius, i.e., it has a maximum for a par- 
ticular value and decreases on both sides of this optimum. Some possible 
relations to biological phenomena were postulated. 

Finally, a short examination was made of the relation of the shape of the 
particle (reaction volume) to the reaction rate. It was shown that for the 
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1-step system the sphere has the slowest rate, but that in 2-step systems 
the sphere can have larger rates than the cylinder. 

On the basis of these considerations it may be concluded that structure 
as studied has significant effects on important biological properties. 


I wish to thank the Committee on Mathematical Biology and, in par- 
ticular, Professor H. D. Landahl and Dr. C. S. Patlak for their help and 
criticism. 

(Note: Appendices IA, IB, IVA, IVB, etc. following are elaborations of 
Sections I, IV, etc. in text.) 


APPENDIX IA 
If \ = 0 in (7), 1 becomes gq. We want to determine the ratio q:/1. 
Since the problem of finding ratios will occur frequently, we shall simplify 
the operation as follows. Let 


Ai +B, _ A+B: 


Wy, sleet = 
See tal W, CoD (78) 
Hence, taking the ratio 
Wi _ (Ai +B,—B2+B,) (at Di— Dit D2) (79) 
W, (Ci+ Dy) (Ai+Bz2) 
or 
UA Gaky yaseet bs Dz) + Ai (D2— Dy) — Be (Dit) (80) 
Ws, (Ci + Dy) (Ai +Bz2) ‘ 
We shall use formula (80) in all future calculations of ratios. 
In this particular case, setting gq: = Wi, 2, = We, we see that 
By =0= D,, A, = Ww Xow, — W_w_2X2 , 
81 
By, = d (w1X0—W-1X-1) ; Ci =w.+w-_ , Dr:=2. oe 


Hence substitution of (81) into (80) results in 


4 os Be hn ee) keel ee 
2% (w2-F wi) [wiXo (a+r) —wyAXp— ww aka § 2) 


Equation (82) is 1 when 
5 as Wi Xo + w_.Xe 
10 = 
Wo +w_y 
it is infinite when 
Rae W,Xo (A+ we) —w_yw_2X> 
W—1Xr 


Since the latter is larger than the former, the conclusions in the text 
follow. 
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APPENDIX IB 
The basic equation to be solved is 
QB 


eae (83) 
with boundary conditions 
dB 
=e DF. ae — k,B (0) ’ 
dB ibe 
= D)= ksB (1) . 
The solution of (83) is given by 
B(x) = + ex. (85) 
The coefficients c, and cz are determined from (84). They are 
be eet ae 
> D (hot ks) + Rakes’ 
(86) 
ae k, A, (D+ kJ) 
‘ D (ho+ ks) + Roksl 
Hence a 
= —— . 87 
BD) — 5G. By) tak i 
(28) follows from (87) directly. 
APPENDIX IVA 
We have 
4 ki 
A+ ae Xi, 
Mie +E 
as (88) 
. Xn= Suti+En, 
with these two conditions 
3 Fu = XT+E; ; A pS 
; (89) 


K= p> (S;+ Xi) +Sn4i- 


Hence 


n 


S,= K=Xi— DSi + Xi) — Sots 
2 * 
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The basic equations describing (88) are 


dX, 
dt 


= kyE,S,;— X;(k-1+ 01) + 0-182 


ae aX — o-1S2 (Eo: — X1) — RoS2 (Eon — Xo) + k-2X2 
(90) 


dSn+i 


dl a= Ce ari T—wSnt+10n + ye, . 
Setting 

a= X;—X;, pal 2 aes 

f= S;—S;, Fas ee rg hs oF 


eliminating all terms in a; or 8; higher than first power, and using (89), 
results in 


d n be be n+1 = 
So ala Xit+ kXi— k,K+ Jo» = Ri Eyo 
eee! fas k_y ia 0-15: | + Be [RX ale Rk, Eyo+ o—1E 9 — ¢-14) 
n+1 4 n = 
— ki >* 8; (Ew — X1) — ki D> a: (Eo — X) , 
3 2 
d = = = 
cP ay (o,+ o-1S2) + Bo (o-1X1 — o—1E 9 — RoEoo + kyXo9) 
+ a2 (keS2+ k-2) (91) 
da; 


“Gp = Baki (Eos — X5) — 05 (0-3Sj41 + b-j+ 03 + 8;S;) 

i + Bj+10~—; (Eo; — X;), 
or = aj;—1 (oj-1 + o-G—-)S;) + B; (o-G_-p Xj — o—(j;—-1) Eoj-1 

— k;Eo; + kj X;) +0; (RS; +k_;). 


It can be seen on inspection of (91) that there is only one pair, a,j, aji, 
which may have opposite signs. This occurs for the coefficient of B. in the 
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first equation and for the coefficient of a; in the second equation. They are 
(Exo — X1) (e..— ky) and 63 -+0-1S2 ) 
respectively. Clearly if k; > o_1, they are of opposite sign. 


APPENDIX IVB 
From (90), if x = 1, 


K=Se+Sit X1; BE, =E\y— Xi, 


225 ky (Eu — X1) (K —S2— X)) 
= Xy (kato) + o-1S: (En — XD; ee 
oo = wy Xy — o= 353 oi — Xi) - 
We set 
a= Xi—-Xi, 
Bo = Sy — Sz 
Hence 
as 1 (Eo — 01 — X1) (K — Bo — S2— 01 — X1) 
— (a, + Xs) (R-1 +01) + 0-1 (B2+S2) (Eo — a1— Xi), 
1s 5, (a, + 1) —0-1(B: +S.) En—a.— Xi). 
Since 


ky (Eo. — X1) (K — S2— Xi) — X,(k-it+ a1) + o1Ss (Eo. — X1) =0 ’ 
re Cat ey (Eo. — X1) =0, 


we find 
oe — ky Eo (Be+ a1) + kX; (ai+ Bo) = kya, (K —S2— Xi) 
mae | (Roi 01) + o—-1B2 (Eo — X,) — o-1S201 y 
os Q101 — o—1B2 (Eo. — X)) ae o—1Sea4. 


We now let a, = me! and az = we"; the above becomes 
Bi [ a A+ k, Eo. — kyXi+ ky es So Smee X;) a [ee ot o—1Sel 
+ ps (Eo: — X1) (ki— o-1) = OF (93) 
—plout o—1S2] tee [— A+ o—1 (En — X1)] =0. 


Equation (44) results from setting the determinant of coefficients of (93) 
equal to zero and expanding. 
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We now calculate 6? — 4ac = G, as defined in the test. It is given by 
G= [o-1(Ei+Se) + hi Gi+51) + 01+ k-1]? (94) 
— 4B, [o-yk-14+ kioi t+ o-1ki (Ax + Si t+ Sz) J , 
where we use the relations 
K=Si+ 52+ X,, 
fu =E,+ Xi. 
Equation (94) can now be expanded and simplified. The resulting equa- 


tion is now written by grouping all terms in the powers of F;. The follow- 
ing results: 


G = Ej (0-1 — hi) 2+ 2B, (0-1 — fx) [o-1S2 — RS + 01 — k=l 
+ (o-1S2+ RS)? + (01+ R-1) ?+2 (01+ b-1) (0-182 + Bi Si) . 
Inspection discloses that G is greater or equal to F, where F is given by 
F = [E,(01— hi) + o-1S2— hiSy + 01 — ka)? 
Since F = 0, G can never be negative. 
APPENDIX IVC 


This system is described by the following set of differential equations 


pee ice tA Ala A.) 07,4 


dl u 
dAjo 
ie mee 3 (Az— A,,.) — k,A,,A 22 TA 45 


GAj ju _, ee 


dt g\4 ii ity) — FA; i Ajar, pi %, +145, 7415 $ (95) 
Aa 
ae aad OE iat as k t 
dt iA;, i414 54, Dau? ee (A iy, j+1 A 4, fret. 
— Fs4y, pex4 gay, jr 
GAanin 


=e A Asa 9) A 


n=l hla” nn n>--nn* 


‘ 
> 
. 
5 
} 
4 
E 
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Following the usual procedure, we set ai; = Ai; — Aij, where Aj; is the 
steady-state value. Hence, ignoring all but terms a;; to the first power, 
(95) becomes 


day, = 
dl =a,,(K—17,—4,,) + 1,%9) 
Goa ya — ag (ty th Ay toy) — kA ge 
A fe loot Ses a Gar ke 11222) 
ee as ee a (r.+k.A ee) 
dt Ie Ie} Ar ae ee) Smad Sea a 7,711 
— O44 518A; p41? 
dajii,j+1 7 (96) 
J J ne as = 
“re Os 541% 4 jt, p41 Oia, gta Tpe1— Tit1, i+ 


ay RA; je at UU preg oats 9+2 


Q 

£ 

3 
a 
) 


dt = n—1,n m—1 ta ie hr oe A) ¥ 


Further, since KAy — n(An — Ay) — ovAy = 0 in the first line in (95), 


— nA 
K 75-0 = a - =Q1. 
11 


Continuing for the second line in (95) gives us 


: A 
hand one che =(p. 


In general, there are 2,10; such that 


Ay it (97) 
ato 
Opie a. 


Vi 


Hence using (96) in (97) we find that the coefficient of \’"-? is given by 


In—1. 


(0 tere Biely 
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and the coefficient of \’"-* is given by 


2n—1 n—1 


ites D) 5 > 0,0;- De- >> es >> RiA;, i414 jar, Pele 


1,j7=1 


Therefore 


2n—1 n—1 n—\ 
Os. oo DIO 2 > hae a REA S414 jai, p41: 
1 


Equation (98) must be less than or equal to zero to satisfy the sufficient 
condition for periodicity. 
In particular, we can substitute (97) into (98). This gives 


me q2. 
a—2a,= (73 tee ry Ae bes) 42> r 
1 7,2 


(98) 


n—1 


— a oo Pray. Crk i+1? 


A 


2A A A A 
= Se cae geet 


Py i x! ine ; 


Clearly, aj — 2a. < 0 if each term in the series obeys this inequality. This 
requires the inequality of (47). 
APPENDIX IVD 


The differential equations describing the system become, in x; and xo, 


d 
ao as (A+ w_-1) xi + Axe, 
(99) 
= AX1— X2 (A+ we) . 
The characteristic equation of (99) is (A + wisry)(A + we — vy) — 2 = 0; 


hence, solving this quadratic equation in 7 results in (49). 


Setting 
X10 = Mit Me X20 = Bi + Be 


from (48), and substituting in (99) gives us 

— bi = — (A+ w-1) mi +ABi 

— Ya (X10— Ha) = — (A+ 4-1) (x10 — Hi) +A (x20— Bi) . 
We can solve (100) for yu: and #;. In particular, equations (50) result. 


\ (100) 
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Further, —(41)/(u2) from (50) is equal to 


la ey ES X10 (¥1 — ¥2) 
He AX20 — X10 (A+ H-1— 1) © 


(101) 


Since, by assumption, y1 > ye, (101) is only greater than 1 if 


AX20 
X10 f 


Aa 


We can now use (49) in the above inequality, remembering that 7 is 
the larger of the two roots. Equation (51) follows immediately. Hence, 
(x20)/(x10) must be less than zero. 


APPENDIX VA 
The basic equations of this system are 
dx ae Rix 2s Roxy 
dt ito x 1+ oxy 


es _kexy | 
pea 


— k3x 


1 


(102) 


We now proceed to find the curves 1 = yi(%) and y2 = ye(x) for which 
dx/dt and dy/dt are zero. It can be seen immediately that for x = 0 = y 
this condition is satisfied. Hence the origin is one of the steady-state 
points of the system. 

Setting each equation of (102) equal to zero and solving for y results in 
the following two equations: 

If dx/dt = 0, then 


ky — kz = k3o4X 
= : 103 
emis aire Oise) hes) os, 
If dy/dt = 0, then 
So eal (104) 
Raoox 


Now consider the case where k, < ks. Therefore the denominator is posi- 
tive and the numerator negative in (103). The resulting curve and its in- 
tersection with the curve (104) is shown in Figure 4. Hence there is only 
one steady state, the origin. This is also stable, as can be seen by examin- 
ing (102). 

If ks < ky < ks + he01/o2, the numerator is positive for small values of 
x, and the denominator is still positive. The resulting graph is shown in 
Figure 5. As can be seen from (104), the origin is no longer stable. 
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If ky > ks + koo1/o2, the numerator is positive for small values of x and 
the denominator is now a quadratic equation in « with two changes in 
sign. Hence it can have two roots, i.e., it can be zero for two values of x. 
The resulting graph is given in Figure 6. As can be seen, four steady states 
are now possible, two of which are stable, two unstable. 


APPENDIX VIA 
We see from (56) that 


ee = Yee (aK +1) =", (105) 
and, from (55) and (58), 


vm eft oof =O] 8 eo HO 


We now approximate each exponential by the first two terms of its Taylor 
expansion. Equation (105) therefore becomes 


Are (« = = 
a Se Oa bie 
Os Ge See eT (106) 
a (Gk -F 1) moj tae = na! 
a A Re =) 
Since 
(Ht) =x =x exp(- 44 a 


(59) follows directly from (106). 


APPENDIX VIB 


The Debye-Hiickel theory for a spherically symmetrical system is 
given by 

dry 

dr? 


= K*ry. 


In the case of the sphere, one boundary condition is given by the require- 
ment that Y remain finite at the origin. This is fulfilled only if 
y ee Wo = Kr ; 


(107) 


We now use (107) to solve for 


dy 


dr}, 


: 
i 
; 
4 
y 
a 
4 
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and set this equal to 41a/e. This becomes 


Ar ao 


e (aK cosh Ka— sinh Ka) ° 


Yo= (108) 


Proceeding as before we find o; and expand the exponentials in a Taylor 
series, ignoring all but the first two terms. Equation (64) then follows. 


APPENDIX VII 


The average number of molecules of type 7 per unit volume of the sphere 
is given by (72). Now 


e—lezw(r)//FT — [1- & = (ORG +5(4 a) ye = 


hence, since 


a OD (n;— Mio) , 
i=1 


P= wy | fren (109) 


x[1- = > (4 ve (r) Jar — mol. 


Or 
Pai 15 ("3 io? <n for (Se —| Ee a Ser] = nok (110) 


Wecan see that the first and last term cancel. Also, from the condition of 
electroneutrality, 


212,021 = 0 ; 


hence the linear term under the integral sign vanishes. The result is 


3 e hy Saar 
— Ni Z ar. 111 
YS wozi f yar (111) 


We integrate (111), using 
y= 73 sinh Kr, 


where <P 
Zz 2 Ale coe 
K=XV (zn) ’ r eee ? 


the result is equation (73). 
Further, let A~ be the anion surface charge of the enzyme, B* the cation 
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surface charge. Let C,; be the ith cation in solution C_; the 7th anion in 
solution. Hence 
K: 
A-~+C4;--> ACH, 


AGES a, 


BtC_,25BC_, 


Bt +C_,23BC_, . 


If A> = the total anionic surface concentration and By = the total cat- 
ionic surface concentration, then 


ee 
14+ K C4: 
: (113) 
B 
BS 
1+20,C-; 
Further, the total surface charge ¢o = B+ — A-. Since 
Car = mpe7W/ET , — C_; = n_;e(W/ET) , 
the surface charge can be written approximately as 
a Brees gees (114) 
DQ;n-; Kini; 
Since 
dy 4no 


dr/, eae 


we can solve for the unknown coefficient Yo. Approximating (a) in (114) 
by the first two terms of the Taylor expansion results in (74). 


APPENDIX VIII 


a(r) =a (1-5), 


C2 (1r) = Ce (1 =% F 


We assume that 


eS —: wo 
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where é, and @ are the concentrations at the center of the sphere of radius 
R. We now equate the total production of each substance to its outflow. 
Hence 


4x f-}K.— Fea ae (1 =) 


(115) 
+ Kats(1 —f)brar =45Ri,D 
and 
arf} KC, (1 —5)- Kb; (1 =) rtd = 4aRCD. 
We use (115) to evaluate C; and C2. We find that 
pee 4 K,R? (K_,R?+ 12D) ; 
144 D?+ 12 DR? (K_,+ K2+ K-2) + K-1K-2R* (116) 


fhe 4K, K2R! 
>= 744 D? + 12 DR® (K-1+ K+ K-2) + K-1K-2R* 


Since the total rate of production of C, is given by ArRC2D,, (76) follows. 
A very similar calculation can be made for a cylinder of radius 4 and 
length zo. We again assume that 


= h 22 
Gties) =6,(1-;) +2 


= h 22 
C, (is) =0.(1-7) ce 


(117) 


Equating the total rate of production of C; and C2 with their outflows 
gives us 


ty hy : : 
an fas f (K+ Ta hed ae VEN 


. (118) 
h 22 An Dh, . ) 
ae ake ae Beale ass D 
x|1 | [1+ | hah t+ C80 
and 
2,/2 “ho = = h We 
an f dz f (KC, - K-01 (1-7) (+2 hdh 
idee ig : (119) 


bee Aa DhiCs 
z 


0 


+ 1C2%0D 7 
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From (118) and (119) C, and C; can be determined. Since 


5 Ah? 2 
Wy = 10s a oly 


0 


(77) follows from the determination of Ce 
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RECEIVED 4-10-54 


Note: It was called to the author’s attention during proofreading that 
there is only one stable point in Figure 6. However, if the following sets 
of numbers are used, there are two stable and two unstable points (the 
graphs are very similar): k, = 21, ke = 20, ks = 1, ka = 6, 0; = 0.6, and 


a ihe 


